Draft version September 28, 2011 

Preprint typeset using I^'T^]X style emulatoapj v. 08/22/09 



O 
(N 

D 
00 

l> 
(N 

6 



> 

o 
o 

On 

vn 
O 



MAPPING THE SHORES OF THE BROWN DWARF DESERT HI: YOUNG MOVING GROUPS 

Evans, T.M.^ ^ Ireland, M.J.^-'"*, Kraus, A.L. Martinache, F. Stewart, P.\ Tuthill, P.G.\ Lacour, S. ^'S 

Carpenter, J.M. * and Hillenbrand, L.A. *^ 

Draft version September 28, 2011 

ABSTRACT 

We present the results of an aperture masking interferometry survey for substellar companions 
around 67 members of the young 8-200 Myr) nearby (~ 5-86 pc) AB Doradus, /3 Pictoris, Hercules- 
Lyra, TW Hya, and Tucana-Horologium stellar associations. Observations were made at near infrared 
wavelengths between 1.2-3.8 /zm using the adaptive optics facilities of the Keck II, VLT UT4, and 
Palomar Hale Telescopes. Typical contrast ratios of ~ 100-200 were achieved at angular separations 
between ~ 40-320 mas, with our survey being 100% complete for companions with masses below 
^^0.25 Mq across this range. We report the discovery of a 0.52 ± 0.09 M© companion to HIP 14807, as 
well as the detections and orbits of previously known stellar companions to HD 16760, HD 113449, and 
HD 160934. We show that the companion to HD 16760 is in a face-on orbit, resulting in an upward 
revision of its mass from M2 sini ~ 14 Mj to M2 = 0.28 ± 0.04 M©. No substellar companions were 
detected around any of our sample members, despite our ability to detect companions with masses 
below 80 Af j for 50 of our targets: of these, our sensitivity extended down to 40 Af j around 30 targets, 
with a subset of 22 subject to the still more stringent limit of 20Afj. A statistical analysis of our 
non-detection of substellar companions allows us to place constraints on their frequency around 0.2- 
1.5 M© stars. In particular, considering companion mass distributions that have been proposed in the 
literature, we obtain an upper limit estimate of ~9-ll% for the frequency of 20-80 A/j companions 
between 3-30 AU at 95% confidence, assuming that their semimajor axes are distributed according to 
dM /da cx in this range. 

Subject headings: stars: binaries: general; stars: low-mass, brown dwarfs; stars: pre-main sequence 



1. INTRODUCTION 



In the past few years, direct imaging surveys have 
begun to build up a picture of the mass and semi- 
major axis distributions of substell ar compa nions at 
separations beyond ~ 20-30 AU (eg. iBiller et al. 2002 
Carson et al. '20061 : iChauvin et all I2 OIOI: 'Kaspcr ct aH 



2007; Lowrancc et al. 



2007: Lafrenic r c et al.l 

Masciadri et al.! ' 120051 : iMetchev fc Hillenbrandl |200£ 



Meanwhile, statistical analyses of radial velocity results 
have tended to focus on objects wit h masses below 
^ lOAfj out to separatio ns of '^3AU l| Gumming et al.l 
I2OO8I : [Howard et al]|2010D . However, given the observa- 
tional biases of radial velocity and direct imaging sur- 
veys, the separation range of '^3-30AU has remained 
relatively unexplored. 

Aperture masking interferometry is a direct detection 
technique that is well suited for detecting substellar com- 
panions with masses of 10 Afj and semimajor axes 
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within '^30AU around young, nearby stars. For in- 
stance, it has been used to conduct surveys for sub- 
stellar compani ons around members of the Uppe r Scor- 
pius (.Kraus et a l. 2008) and Taurus- Auriga (.Kraus et al.l 
1201 Ih associations, as well as to measure the dynam- 
ical mass of t he brown dwarf companion to GJ 802b 
([Ireland et al.l UOOS). show that CoKu Tau/4 is a binary 
syst em rather than a transitional disk (llrela nd fc KrausI 
1200 8). and place limits on pos sible cornpanio ns existing 
within 10 AU of HR 8799 (|Hinklev et al.ll2011[ ). Recently, 
the technique has also produced the first direct detection 
of a young exoplanet still undergoing formation within 
the transitional disk of LkGal5 (Kraus & Ireland, sub- 
mitted) and a similar detection of an objec t within the 
gap of the T Cha disk (iHu elamo eraLll2011[ ). 

This paper presents the results of an aperture mask- 
ing survey of 67 members of the AB Doradus (AB Dor) , 
/3 Pictoris (/3 Pic), Hercules-Lyra (Her-Lyr), Tucana- 
Horologium (Tuc-Hor), and TW Hya (TWA) moving 
groups. At least 49 of our targets have been observed 
previously as part of deep imaging surveys, but these 
observations have typically been sensitive to different or- 
bital separations than those that are probed here. We 
chose our targets based on their youth (8-200 Myr) and 
proximity (5-86 pc). The former ensures that any sub- 
stellar companions are still glowing relatively brightly 
at infrared wavelengths following their recent formation, 
while the latter allows smaller absolute separations to be 
explored for a given angular separation. 

The paper is organized as follows. In Section [21 we 
provide a brief overview of the aperture masking tech- 
nique. In Section [31 we describe our survey sample. In 
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Section 2] we summarize the observations that were made 
and how the data were reduced. In Section [5l we explain 
how we searched for companions to the target stars in 
the reduced data and how we derived the survey detec- 
tion hmits. In Section [6] we report our resuhs, including 
the detections and orbits for stellar companions around 
HIP 14807, HD 16760, HD 113449, and HD 160934. How- 
ever, no substellar companions were detected, and in Sec- 
tion [7] we present a statistical analysis of this null result 
before concluding in Section [H 

2. APERTURE MASKING 

The aperture masking technique works by placing an 
opaque, p erforated mask at or near the p upil plane of a 
tele scope (iFizeaTjl lSGS. 'Michelson|l l891at more recently, 
see iTuthill et arm OOO. 2006, 20101). This converts the 
single aperture into a multi-element interferometer. Each 
pair of holes in the mask acts as an interferometric base- 
line, resulting in an interferogram being projected onto 
the image plane. 

The complex visibility V (jMichelsonI Il891bf ) of the 
source brightness distribution S is sampled by taking 
the 2D Fourier transform of the measured interfero- 
gram /. This follows from the Van Cittert-Zernike Theo- 
rem, which states that the normalized complex visibility 
is equal to the Fourier transform of the brightness distri- 
bution: 



V = 



So 



(1) 



where the tilde denotes the Fourier transform and So is 
the total source flux. Since the detected image is the con- 
volution of the instrumental point spread function (PSF) 
and the source brightness distribution, this leads to: 



SnP 



(2) 



where P denotes the Fourier transform of the PSF. In 
practice, the PSF is measured by observing an unresolved 
calibrator star, i.e. a point source, which has unit com- 
plex visibility V — 1. 

In this study, we used non-redundant aperture masks, 
with each baseline pair corresponding to a unique point 
on the spatial frequency plane. We used masks with 7, 
9, and 18 holes, giving 21, 36, and 153 independent base- 
lines, respectively. Hole diameters and transmission frac- 
tions are provided in Table [TJ The subaperture configu- 
rations on the masks were designed to provide a uniform 
and isotropic sampling of the complex visibility function, 
with the specific mask chosen to observe a given target 
depending on the target's brightness and the expected 
sources of systematic error. For example, although the 
18 hole masks had slightly longer baselines than the 7 
or 9 hole masks, they had lower total throughput with a 
broader PSF. This meant that they could only be used 
with narrow band filters, which restricted their use to 
brighter targets. 

To identify faint companions around our targets, 
we used a quantity derived from th e complex visi- 
bility known as the closure phase Q ([Jenniso i [19581: 
[Baldwin et al. I [19861) . The closure phase is obtained by 
adding the complex visibility phases around a closure 
triangle of subapertures. Explicitly, if we denote the 



measured complex visibility phase between the zth and 
jth subapertures as (fij, the intrinsic complex visibility 
phase as 4>ij, and a phase error due to atmospheric and 
instrumental effects across the zth aperture as ?7i, then 
we have: 



fij = (fiij +r]i - Vj 
^ij = (l>jk + - Vk 



(3) 



Importantly, the diameter of the mask holes are chosen 
to ensure that the wavefront phase variations across each 
subaperture are approximately constant so that they 
can be neglected. Combining aperture masking with 
adaptive optics allows subaperture diameters that are 
larger than the atmospheric Fried parameter and expo- 
sure times that are longer than the atmospheric coher- 
ence time to be used, providing a greater through-put of 
photons. It follows that the r]i terms cancel out when we 
take the closure phase sum: 



6ijfe = 4>ij + 4>]k + 't'ki 



(4) 



where Qijk is the closure phase of the triangle ijk. 

The independence of the closure phase quantity from 
major sources of systematic error allows us to achieve the 
full interferometric resolution according to the Michelson 
criterion, which is equal to X/2B, where A is the ob- 
serving wavelength and B is the longest baseline on our 
mask. This is the smallest angular separation for which 
two point sources would be fully resolved. Given that 
the longest baseline of the masks used in this study span 
nearly the entire telescope aperture, this corresponds to 
angular scales of roughly half the single-aperture diffrac- 
tion limit. 

3. SURVEY SAMPLE 

In 2007, our group initiated a search for close, faint 
companions around young, nearby stars using the aper- 
ture masking facilities installed on the 5.1m Hale tele- 
scope at Palomar Observatory in California. In subse- 
quent years (2007-2011) the survey was extended and 
made use of similar facilities installed on the 10 m Keck II 
telescope at Keck Observatory in Hawaii and the 8.2 m 
VLT UT4 telescope at the VLT Observatory in Chile. 

Our final survey sample consist ed of 67 p r opose d 
memb er s of the AB Dor ( |Z uckemian et al.[ I2004D. 
Pi c (|Zuckerma.n et al.l l2001aD. Her-Ly i~( FuhrmannI 
[200l. Tuc-Hor (iZuckerman et all l2001b[ ). and TWA 
([Kastner et al.lll997D moving groups. A concise summary 
of the sample is provided in Table [2| while the full list is 
given in Table [31 Figure [T] shows the sample members 
binned according to spectral type and masses. 

In selecting our targets, we noted that many of the 
moving group members have already been identified as 
binary systems. The presence of a binary companion 
within ^ 1" of a target star reduces the ability of aper- 
ture masking to detect additional companions, because 
the interferograms will overlap. Also, similar brightness 
companions at separations of ^1.5-3" can prevent the 
adaptive optics system from achieving a stable lock on 
the target. For these reasons, we chose not to include any 
targets in our sample that were known to be affected by 
such issues. 
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Fig. 1. — The 67 survey targets binned according to the spectral 
types (top panel) and masses (bottom panel) listed in Table [S] 

We also emphasize the difficuhy of assigning moving 
group membership to individual stars. Consequently, it 
is possible that not all objects in our sample are neces- 
sarily young. In particular, the moving group member- 
ships of nine of our targets (HD 89744, HD 92945, GJ466, 
EKDra, HIP 30030, TWA-21, TWA-6, TWA-14, TWA- 
23) were eit her unable to be con firmed or else ruled to be 
unlikely by iTorres et all (pOOl ) using a dynamical con- 
vergence analysis. Furthermore, the existence of Her-Lyr 
as a genuine moving group is not yet as well-established 
as the others. To investigate how sensitive our statisti- 
cal analysis presented in Section 17.21 is to the uncertain 
membership of these targets, we repeated the calcula- 
tions separately with them included and then removed 
from the sample. 

4. OBSERVATIONS AND DATA REDUCTION 

We observed our program objects over the course of 
twelve observing runs using the facility adaptive optics 
imagers at Palomar (PHARO), Keck (NIRC2), and VLT 
(CONICA) between April 2007 and April 2011. Each 
camera has aperture masks installed at (Palomar, VLT) 
or near (Keck) the pupil stop wheels. The central wave- 
lengths and bandpass widths for each filter used are listed 



in Tableland details of our observations are summarized 
in Table [5l Observing conditions varied widely, but we 
attempted to match the observations to the appropriate 
conditions. Most of our brighter targets were observed 
through clouds or marginal seeing as they were the only 
ones we could lock the AO system on, while our fainter 
targets were typically observed under better conditions. 

Our observing stra tegy has been described previously 
in lKraus et al.l ()2008() . Each observation consisted of 1-3 
target-calibrator pairs, usually with ~ 10-20 frames per 
block. We tried to choose calibrators with optical and 
near-infrared brightnesses that were similar to those of 
the target, rather than calibrators that were necessar- 
ily brighter. This was done due to concerns about the 
magnitude-dependence of non-common path errors in the 
adaptive optics system. For targets of brightness R ^7, 
calibrat ors were chose n from the stable radial velocity 
stars of INidever et al.l ([2002) . For fainter stars, we could 
not explicitly choose calibrators that had been vetted for 
close binaries, so we simply chose nearby 2MASS sources 
with similar colors and brightnesses. In all cases, we tried 
to select calibrators that appeared to be single and un- 
blended in the 2MASS images, as well as close to the 
target on the sky (< lOdeg for the Nidever et al sources 
and <3deg for the 2MASS sources). In addition to re- 
ducing overhead times, using nearby calibrators helped 
to minimize residual wavefront errors introduced by long 
telescope slews. 

Data reduction was performed using our group's 
custom-written IDL pipeline (for fu rther details, see 
Llovd et al.l 120061 llreland et ahl 120081 and iKraus et al.) 
20081) . Complex visibilities were extracted by Fourier- 
inverting the cleaned data cubes and sampling the uv- 
plane at points corresponding to the mask baselines. 
Calibration was performed by subtracting the calibra- 
tor complex visibility phases from the complex visibility 
phases of the science targets. 

5. BINARY MODEL FITTING 

We used the same method as IKraus et"all ()2008l[2011| ) 

to search for companions to our targets over the sepa- 
ration range 20-320 mas. We only used closure phases 
in our binary model fitting, discarding the visibility am- 
plitudes as these are more affected by systematic errors. 
The parameters we fit for were the angular separation 
p between the primary and companion, the position an- 
gle 9 of the companion, and the brightness contrast ra- 
tio C = fp/ fc, where fp and fc were the fluxes of the 
primary and companion, respectively. Fitting was per- 
formed by initially fixing a high contrast ratio of C = 250 
and generating the corresponding model closure phases 
for each point on a grid of angular separations span- 
ning 20 < p < 320 mas and position angles spanning 
< 9 < 360 deg. The point on the p-9 grid giving 
the lowest for the measured closure phase values was 
then taken to be the starting point for a steepest-descent 
search in which all three model parameters (C, p, 9) were 
allowed to vary. The initial grid search ensured that the 
final minimum reached corresponded to the global mini- 
mum. 

The binary fit was considered to be bona fide if it 
passed a 9 9.5% detecti o n criterion, which has been ex- 
plained in IKraus et"all ()2008l 120111 ). This was done by 
generating 10 000 artificial closure phase data sets with 
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Fourier plane sampling that was identical to that of the 
measured data. Each artificial closure phase was ran- 
domly sampled from a Gaussian distribution with a mean 
of zero, corresponding to an unresolved point source, and 
the same variance as the corresponding measured value. 
A best-fit companion contrast C was then obtained for 
each set of artificial closure phases using minimization 
at each point on the p-9 grid. Once again, by searching 
over the entire grid we ensured that the global minimum 
was identified. A 99.9% detection threshold was then 
defined separately for five contiguous annuli (20-40 mas, 
40-80 mas, 80-160 mas, 160-240 mas, 240-320 mas), cor- 
responding to the 0.1th percentile of the best-fit con- 
trasts obtained for the artificial data sets within that 
annulus. In other words, if the target was a point source 
instead of a binary, there was only a 0.1% chance that 
the measured closure phases would give a best-fit con- 
trast lower than the threshold value in the annulus cor- 
responding to the best-fit separation. This corresponds 
to a 5 X 0.1% = 0.5% false alarm probability across the 
full 20-320 mas range. Therefore, if the best- fit model 
satisfied this condition, the detection was considered to 
be real at 99.5% confidence. 

It was important to ensure that any high probabil- 
ity (>99.5%) detections were not caused by companions 
around one of the calibrators rather than around the sci- 
ence target. A small number of such false alarms (~5) 
did occur during the course of our analysis. Such cases 
were usually quite straightforward to identify by system- 
atically repeating the calibration and binary fitting, ex- 
cluding one calibrator at a time. Given that the calibra- 
tors did not have known ages, but were likely to be Gyr 
old, any companions detected around them were almost 
certainly not substellar, and so they were not considered 
further. 

6. RESULTS 

Using the method described in Section [SJ we identi- 
fied stellar companions to four of our AB Dor targets 
(HIP 14807, HD 16760, HD 113449, HD 160934) and re- 
port our best-fit binary solutions in Table [6] Of these, 
the companion to HIP 14807 is a new discovery while 
the companions to HD 16760, HD 113449, and HD 160934 
are the same as those di scovered independently us- 
ing radial velocity ( Bouchv et al.ll2009t I Sato et al.l l2"00l 
iCusano et ahl [200l i2010i: iGalvez et al.l 120061 ). We de- 
scribe the detected companions in Sections 16.1146.41 a nd 
present our full survey detection limits in Section [6.51 

6.1. HIP 14807 

A companion was clearly detected in our Keck observa- 
tions of HIP 14807 on 2009 November 21 (MJD 55156.2) 
at an angular separation oi p — 28.74 ± 0.19 mas with a 
contrast ratio of C = 3.00 ± 0.06 in the CO fiher. As- 
suming a system age of 1 10 zfc 40Myr, i nterpo lation of 
the NextGen isochrones of iBarafFe et al.l ()1998[ ) gives an 
estimated companion mass of 0.52 ± 0.09 Mq, which in- 
cludes the uncertainty in the age and distance, as well as 
the uncertainty in the fitted contrast. 

The companion was also detected at high confidence 
in the Palomar data from 2007 November 29 (MJD 
54433.1), with a fitted contrast ratio of C = 10.15 ±3.71. 
However, this error bar is neither symmetric nor real- 
istic, as there is a strong degeneracy between contrast 
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Fig. 2. — The contrast /separation degeneracy for the fit to the 
HIP 14807 data taken in K band on 2007 Nov 29 (MJD 54433.1), 
showing the default fit from our pipeline (triangle) and the fit after 
fixing the contrast (star). The third axis of position angle has been 
integrated over to give a two-dimensional likelihood plot. Contours 
are nominally at 90, 99, 99.9, and 99.99% confidence. 

and separation for small separations in aperture mask- 
ing data se ts. This is illustra. ted in Figu re [2] (see also 
Figur e 7 in iPravdo et al]l2006L Table 2 in llreland et aP 
mm. and Fi gure 2 in iHuelamo et"alll2011l) . A fuller dis- 
cussion of this degenera cy is provided in Section 2.1 of 
iMartinache et all ( 20090 . In cases such as these, quick 
data sets were taken with only one or two calibration 
observations. As a result, the quoted error bars are not 
necessarily accurate at the few tens of a percent level, be- 
cause the dispersion between calibrators is used to esti- 
mate the errors in the closure phases. Despite this, global 
orbital fitting to multiple aperture masking data sets has 
been performed successfully by using a single contrast 
for all epochs, with the resulting astrometric fits being 
consistent with those obtained using othe r techniques, 
and having reduced of order unity (eg. IDupuv et al.l 
[2009b). 

For these reasons, we redid the binary fitting to the 
MJD 54433.1 data using a prior on the contrast deter- 
mined from the other well-constrained fit to the MJD 
55156.2 data. The results of this revised fit are given in 
Table [71 

6.2. HD 16760 

HD 16760 is unusual because it shows signs of be- 
ing both young and old. Its youth is implied by 
its high lithium abundance, as well as its physi- 
cal association with the active star HIP 12635 and 
a common proper motion with the AB Dor group 
([Zuckerman et al.l [20041: [Torres et al1[2008D . Its old age 
is implied by its low vsvai value (2.8 ± 1.0 kms~^, 
IBouchv et aL[[2009t 0.5 ± 0.5 kms-^ [SatoeFall [20091) 
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and its low Ca H fc K activity index (logi?^^^- = 
-4.93. ISato et al l [20091 logi?^^ = -5.0 ± 0.1, 
iBouchv et al.ll2009D . which is consistent with field dwarfs 
(logi?^^ = -4.99 ± 0.07. iMamaiek fc Hillenb rand 200|) 
and inconsistent with high probability members of the 
625 Myr Hyades cluster (loe R'„^ = -4.47 ± 0.09, 
IMamaiek fc Hillenbrandl 120081) and other youn g stars 
(see Tables 5-8 of IMamaiek fc Hillenbrandl [20081 ). How- 
ever, we note that this system is not the only example 
of a binary pair showing contradictory age indicators: 
when examining the activity consistency o f known bi- 
nary pairs, IMamaiek fc Hillenbrandl (|2008[ ) identified a 
similar case of an inactive primary with an active com- 
panion (HD 137763 A/B). 

Previous radial velocity measurements h ave shown 
that HD 16760 possesses a close companion (jSato et al.l 
120091 : IBouchv et all [2009 1). for which Sato and cowork- 
ers derived a minimum mass M2 sini value of 13.13 ± 
0.56 Mj, while Bouchy and coworkers obtained a sim- 
ilar value of 14.3 ± 0.9 Mj. We clearly detected this 
companion in our Keck data from 2008 December 23 
(MJD 54823.2), 2009 August 6 (MJD 55049.6) and 2009 
November 20 (MJD 55155.3) (Table El- Taking the 
weighted mean of the well-constrained isochrone mass 
estimates, we obtain a mass of M2 = 0.28 ± 0.04 M© 
for the companion, which includes the uncertainty in the 
age, distance and fitted contrasts. This places it well 
within the stellar mass range. 

Meanwhile, due to degeneracy between contrast and 
separation (see Section 16. 1|) . combined with mediocre 
data quality, the separation derived from the MJD 
55155.3 K band data was inconsistent with the separa- 
tion derived from the J and H band data taken on the 
same night. For this reason, we obtained a further obser- 
vation the following night with the CO filter, which has 
a very similar bandpass to the Kcont filter (see Table [4} 
we had intended to use the Kcont filter, but there was 
a mix-up in the filter selections). The binary parame- 
ters derived from this follow-up observation are in close 
agreement with the values obtained from the J and H 
band observations. 

We also repeated the fit to the degenerate K band data 
with a prior on the contrast obtained by combining the 
fitted contrasts to the other K band epochs. The system 
properties derived from this re- analysis are reported in 
Table [3 and agree with the values obtained for the J and 
H band data sets. We note, however, that the calibra- 
tion error added in quadrature to obtain a reduced 
of unity for this fit was 1.8 degrees. This is unusually 
large and suggests that the quoted errors for the inferred 
parameters are likely to be underestimated somewhat. 

Using our multi-epoch data, we were able to derive an 
orbital solution for the companion. To do this, we fixed 
the values for the time of periastron Tq, orbital period P, 
orbital eccentricity e, and argument o f periastron uj pub - 
hshed for the radial velocity orbit from lSato et al.l ()2009[ ). 
We were not able to fit for the orbital inclination i using 
our aperture masking astrometry data because we only 
measure the axis ratio of the visual orbit, and this varies 
with the cosine of the inclination. Hence, we are not sen- 
sitive to small changes in i when i is near zero, as is the 
case here. Instead, we combined the model-dependent 
mass estimate obtained from the aperture masking re- 
sults with the value for M2 sin i derived from the radial 



o 

u 

0} 

a 




30 20 10 -10 -20 -30 

RA Offset (mas) 

Fig. 3. — Plotted orbital solution for the companion to HD 16760. 
The black crosses show the aperture masking astrometry points 
with associated uncertainties and the red triangles mark the cor- 
responding epochs of the orbital solution. 

velocity results to calculate i — 2.6± 0.5 deg. Then with 
these parameters held fixed, we inferred values for the 
longitude of the ascending node f2 and semimajor axis a 
by fitting to the aperture masking astrometry listed in 
Tables [6] and [3 The final orbital solution is reported in 
Table [SI and plotted in Figure [3] 

Lastly, we note that the rotational velocity of the pri- 
mary is revised upwards from vsini ^ 0.5-4 kms^^ to 
V ^20-25 kms~^, a value that is more in line with other 
members of AB Dor. However, the low Ca H fc K emis- 
sion of HD 16760 remains unexplained. The only reason 
we might expect to see an inclination dependence in the 
strength of log i?^^ is if the Ca H & K emission varies 
with latitude on a star, such that polar areas show lit- 
tle emission compared to equatorial regions. We are not 
aware of any model that would predict this effect. 

6.3. HD 113449 

We detected a companion around HD 113449 in six of 
our data sets taken at Palomar and Keck between 2007 
Apr 6 and 2010 Apr 25 (Table [6|). Four of these data 
sets (MJD 54196.3, MJD 54634.2, MJD 54821.7, MJD 
553311.4) were well-constrained, and taken together im- 
ply a companion mass of 0.51 ± 0.01 M© b ased on the 
NextGen isochrones of iBaraff e et al.l (|1998l ) , including 
the uncertainty in the age, distance and fitted contrasts. 
The other two data sets, however, gave fits that were 
degenerate in contrast and separation, as has been de- 
scribed above. For the first of these (MJD 54197.6), we 
repeated the analysis with a prior on the separation taken 
from the well-constrained fit to the previous night's data 
(see footnote in Table [6]). For the second case (MJD 
54252.1), the analysis was repeated with a prior on the 
contrast taken from the well-constrained solutions ob- 
tained for the four other H band data sets (see Table [7]) . 

The companion we r e port here was first announced 
by 'Cusan o et al.l (|2009l . I2010( ) subsequent to the com- 
mencement of our survey. Using radial velocity measure- 
ments, those authors obtained a value of F{Mi, M2, i) = 
0.0467 ± 0.0006 Mq for the spectroscopic mass func- 
tion and estimated a secondary-to-primary mass ratio 
of g = 0.57 ± 0.05. In addition, using astrometry mea- 
surements made with the VLT-I they obtained a value of 
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i = 57 ± 3° for the inclination, 17 = 124 ± 4° for the lon- 
gitude of the ascending node and a — 0.750 ± 0.030 AU 
for the semimajor axis. 

We computed an orbital solution for the companion 
using our aperture masking astrometry (Tables IHl and [7]), 
allowing i, a, and f2 to vary as free parameters in our fit- 
ting, while holding P, Tq, e and w fixed at the values de- 
termined Cusano and coworkers. However, we found that 
we could not obtain a reasonable value with the period 
of P = 215.9 days reported by those authors. Instead, 
an acceptable fit was made when we allowed the period 
to be a free parameter, obtaining P — 216.9 days. The 
best-fit parameters are reported in Table [5] and the cor- 
responding orbit is plotted in Figure |4l In particular, our 
fitted value of £7 = 202.0 ± 1.6° d oes not agree with th e 
value of f2 = 124 ± 4° reported in iCusano et al.l (|2010( ). 
but as the details of those VLT-I observations are not 
given, we cannot make a further comparison. Lastly, the 
dynamical mass of the system (Mtot = 1-10 ± 0.09 Mq) 
appears to be underestimated by ~ 2a when compared to 
the isochrone-determined masses (Mi = 0.84 ± 0.08 Mq 
and M2 = 0.51 ± 0.01 Mq). As the orbital period is ~1 
year and the astrometric semimajor axis is comparable to 
the parallax, examining this discrepancy in more detail 
would require refitting to the raw HIPPARCOS data. 

6.4. HD 160934 

We detected a companion around HD 160934 in our 
Palomar data taken on 2008 June 23 (MJD 54640.3) and 
Keck data taken on 2010 April 26 (MJD 55312.6) and 
2011 April 23 (MJD 55674.6). The binary solutions are 
all in excellent agreement (Table [5]). We obtain a value 
of 0.54 ± 0.01 Mq for the companion mass by combining 
the estimates from each epoch. 

Thi s companion was first reported by iGalvez et al.l 
(|2006| 1. who identified HD 160934 as a spectroscopic bi- 
nary with an estimated period of ^^17.1 years and an 
eccentricity of e'~0.8. However, these were preliminary 
values based on limited phase sampling, and a period 
of approximately half this is also consistent with the 
da ta. In fact, this shorter period is the one preferred 
bv lGriffin fc Filiz AkI (j2010( ). who repeated the fit to the 
same data with a small number of more modern radial 
velocity measurements. 

Furthermore, in addition to our values presented in 
Table [6l rel ative astrom etry measurements hav e been 

tblished bv lHormuth e t al. (20Q3) and [LafrcnicreleLalJ 
)07). Using the combined data set, which is summa- 
rized in Table [SI we performed a least-squares orbital fit 
and report the results in Table [51 The solution is plotted 
in Figure [5] In order to achieve a reduced of 1.0, we 
had to add an extra position angle error of 0.4 degrees to 
all data, which may indicate a small position angle cal- 
ibration mismatch between the three instruments used 
in this fit. Of these parameters, only Tq has an uncer- 
tainty that would be significantly changed by the addi- 
tion of radial velocity data, which have not been made 
available to us because at least one new paper including 
those data is in preparation (Monies, private communi- 
cation). However, we can combine th e semiamplitude of 
the ra dial velocity curve published in iGriffin fc Filiz A^ 
(|2010D (ii:i=7.39±0.22kms-i) with our o rbital fit and 
the H IPPARCOS parallax of 30.2±2mas ([van LeeuwenI 
|2007() to obtain a mass of O.48±O.O6M0 for the compan- 
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Fig. 4. — The same as Figure (3] showing the orbital solution for 
the companion to HD 113449. 



100 




-100 f 1 

100 -100 -200 -300 

RA Offset (mas) 

Fig. 5. — The same as Figures [3] and [J] showing the orbital so- 
lution for the companion to HD 160934. Nominal error bars that 
are too small to be accurately represented are instead contained 
within filled squares. 

ion. This value is consistent with the one derived above 
using isochrones, at the level of the uncertainties. 

Although the binary orbit is not taken into account 
in computing the HIPPARCOS parallax, the period is 
several times longer than the length of the HIPPARCOS 
mission and the system was near apastron at the time of 
observations, so we do not expect the orbital photocenter 
motion to have a significant effect on the measured par- 
allax. As the parallax uncertainty dominates our mass 
uncertainty, we have repeated the orbital calculation at 
several fixed parallax values as gi ven in Table [T Ol Ac- 
cording to the NextGen models of lBaraffe et al.l (|1998). 
plausible ages for the companion range from '^50Myr 
through to the zero-age main sequence. Therefore, the 
dynamical mass does not allow us to place a strong con- 
straint on the system age, but the lower range of allowed 
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values is compatible with the age of AB Dor. 

6.5. Survey Detection Limits 

We list our detection limits in Table [Til correspond- 
ing to the 99.9% threshold values for each of the sep- 
aration annuli, as defined in Section [5] These were 
translated into upper limits for companion masses by 
first converting the contrast ratios into absolute com- 
panion magnitudes using the distances listed in Table [3l 
Then combining these intrinsic luminosities with the as- 
sumed ages listed in Table [21 we determined the corre- 
sponding companion mass by interpolating an appropri- 
ate set of isochrones: specifically, the DUSTY isochrones 
(|Chabrier et al.l[2000l ) for objects wi th 1400 K < T^ff < 
2800 K and the NextGen isochrones (IBaraffe et al.ll 19981) 
for objects with Tcff' ^ 2800 K. For the four targets with 
detected companions (HIP 14807, HD 16760, HD 113449, 
HD 160934) we quote the limits obtained for the residual 
closure phases. 

It should be emphasized that the mass limits quoted in 
Table fTTl inherit the s ystematic errors of t he models used 
to compute t h em (e g. lBaraflFe et al]|2002D . For instance, 
iMarlev et al.l (j2007D have shown that the predicted lu- 
minosities are highly dependent on the treatment of ini- 
tial conditions, with "cold start" models generating lu- 
minosities that can be orders of magnitudes fainter than 
those obtained by the "hot start" DUSTY models over 
Gyr time scales. However, objects in the mass range 
that our survey is sensitive to would most likely have 
formed by the gravitational collapse of instabilities in 
the stellar disk, a process that is more akin to the hot 
start scenario. Indeed, recent observational evidence ap- 
pears to fav or the hot start mod els down to masses of 
~10Mj (eg. ILagrange et al.l [20100 or even suggest that 
they could even ov erpredict the luminosity of such ob- 
jects (jPupuv et al.l '2009a. 2010). In the latter case, the 
values quoted in Table [TT] would be conservative esti- 
mates. 

With these considerations in mind. Figure [Hi shows the 
detection limits plotted in terms of equivalent companion 
mass as a function of angular separation. Due to the het- 
erogeneous nature of our observations, which were made 
using different instruments with different filters, we have 
divided the targets into three groups in these plots. The 
top panel shows the detection limits for our older AB 
Dor (~110Myr) and Her-Lyr (--200 Myr) targets, the 
middle panel shows the detection limits for our younger 
/3 Pic ('-12Myr) and Tuc-Hor (~30Myr) targets, and 
the bottom panel shows the detection limits for the TWA 
(^8Myr) targets. The TWA targets have been plotted 
on their own because all but three of them were observed 
during the same observing run at VLT using the L' filter 
with a 7-hole mask. 

7. SUBSTELLAR COMPANION FREQUENCIES 

We have used our detection limits listed in Ta- 
ble [11] to constrain the frequency of ~ 20-80 Mj 
companions in ~3-30AU orbits around 0.2-1.5 
stars. To do this, we employed the same method- 
ology as [C arson et al. (2006), Lafrcniere et al. (2007), 
[Nielsen et al.i (i200a). ,Metchev fc Hillenbrand (20091 
and iChauvin et aL 

nioiS)- We present a brief outline 
of the approach here, but the reader may consult those 
works for further details. 
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Fig. 6. — Detection limits as a function of angular separation 
for the AB Dor and Her-Lyr targets (top), the /9 Pic and Tuc-Hor 
targets (middle), and the TWA targets (bottom). In each plot, 
solid black lines show the median detection limits within the 20- 
40, 40-80, 80-160, 160-240, and 240-320 mas annuli, dashed black 
lines show the 25th and 75th percentiles, and solid gray lines show 
the individual detection limits. 
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7.1. Mathematical Framework 

Firstly, if we denote the outcome of our survey of Ng 
stars as the set {dj}, where dj is equal to zero if no 
companion was detected around the jth star or equal 
to one if a companion was detected, then the likelihood 
that the fraction of stars with companions is equal to / 
is given by: 

Ci{d,}\f) P{f) 



^(/IK}) = 7T 



J,Ci{d,}\f) P{f)df 



(5) 



where C {{dj}\f ) is the likelihood of our data and P[f ) 
is the prior probability that the underlying companion 
frequency is equal to /. We adopt an ignorant prior of 
P{f) = 1- 

The fact that we did not detect any 20-80 Mj com- 
panions allows us to place an upper limit /,i on their 
frequency by integrating Equation [Sj such that: 

a = —, (6) 

!^C[{d,}\f)df 

where a is a fraction giving the confidence of our limit 
(eg. a = 0.95 corresponds to a confidence of 95%). 

Using Poisson statistics, it can be shown that a null 
result implies: 



^(K}i/)=n^ 



(7) 



where pj is the probability that a substellar companion 
would have been detected around the jth star if there 
had been one present. 

7.2. Monte Carlo Analysis 

The next task is to determine values for each of the pj 
terms, and we did this using a Monte Carlo (MC) ap- 
proach. For each target star in our sample, we generated 
10 000 hypothetical companions, each with a mass M2 
and angular separation p. The companion masses were 
either obtained by randomly sampling from an appropri- 
ate distribution (see Section [731 below ) or else they were 
set to a fixed value (see Section 17.61 below) . To obtain 
the angular separations, we had to properly take into ac- 
count the companion orbital eccentricities, phases, and 
ori entations. We did this using the approach described 
bv iBrandeker et all (120061 ) in their Appendices 1 and 2. 
As with the companion masses, this required either ran- 
domly sampling these properties from appropriate distri- 
butions (see Sections 17.41 and 17.51 bel ow) or else setting 
them to fixed values (see Section [7?6l below) . 

Having generated 10 000 hypothetical companions with 
masses and angular separations for each of the targets in 
our sample, we then consulted the detection limits in Ta- 
bled!] Companions with masses that fell above the min- 
imum detectable mass in the corresponding separation 
annulus for their target star were counted as detections. 
The Pj value for each target was thus given by the num- 
ber Xj of such detections divided by the total number of 
hypothetical companions generated, i.e. pj = a;j/10 000. 
Equipped with the pj values, we were then able to cal- 
culate an estimate for the companion frequency upper 
limit fu at some level of confidence a by integrating 
Equation [6] 



7.3. Mass Distributions 

The distribution of substellar companion masses in the 
separation range ~ 3-30 AU is not yet constrained by ob- 
servations. To accomodate this uncertainty, we have re- 
peated our MC analysis separately for three different as- 
sumed forms for the mass distribution. For the first of 
these, we extrapolated to 20-80 Mj the power law dis- 
tribution that has been uncovered by the Keck radial 
velocity survey for compani ons with masses -M2 < 10 Afj 
and periods P < 2000 days (jCumming et al.|[2008[ ). given 
by: 



dMo 



1.31 



(8) 



where dAf is the number of objects with masses in the 
interval [M2,M2 + dM]. 

The second distribution that we used was the univer- 
sal mass function proposed by iMetchev fc Hi llenbrandl 
(j2009) for companions to solar mass stars, suggested by 
those authors for companion masses between 0.01 Mq 
and 1.0 Mq and semimajor axes between OAU and 
1590 AU. It is given by: 



dAf 
dq 



cx q 



-0.39 



(9) 



where q is the secondary-to-primary mass ratio. 

The last distribution that w e used was a log-n ormal 
paramaterization proposed by iKraus et al.l (l2008f ). de- 



rived using an ad hoc physical model of binary formation. 
It is given by: 



dAf 1 
—— cx - exp 
dq q 



logio Q 



(10) 



and we used the authors' proposed value of cr = 0.428. 

7.4. Eccentricity Distributions 

As with masses, the distribution of substellar compan- 
ion orbital eccentricities in the semimajor axis range 3- 
30 AU is not yet constrained by observations. We chose 
to draw companion eccentricities from a distribution of 
the form: 



dAf 

de 



cx 2e 



(11) 



which, as noted in Appendix 2 of lBrandeker et all ()2006( ) , 
is an approximation that has been derived from physical 
considerations. 

However, to test how sensitive our results were to the 
distribution of companion eccentricities, we repeated all 
of our MC analyses for two limiting cases: (1) fixing the 
orbital eccentricity of all hypothetical companions to e = 
0.9; (2) fixing all hypothetical companion eccentricities 
to e = 0. 

7.5. Semimajor Axis Distributions 

We drew substellar companion semimajor axes from an 
inverse power law of the form: 

dAf 



da 



cx a 



(12) 



over the separation ra nge 3-30 AU. Th i s dist ribution 
was also used bv .Metchev fc Hillenbrandl (|2009[) for a > 
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Fig. 7. — Curves showing the pj values for the case of semimajor axes distributed according to dAf/da oc (Equation I12II and 

eccentricities distributed according to dM jde oc 2e (Equation [11} • The left panel shows the pj values obtained for fixed companion masses 
of 20 Mj (solid lines), 40 Mj (dash-dot lines), and 60 Mj (dashed lines), while the right panel shows the same values obtained for three 
assumed companion mass distributions: a mass power law given by Equation [8] (soli d li nes), a mass ratio power law given by Equation [9] 
(dash-dot lines), and a log- normal mass ratio parameterization given by Equation 1101 (dashed lines). In both panels, the black curves 
indicate the values obtained using the aperture masking limits only, while the red curves show the equivalent values obtained when the 
imaging limits are also included (see Section 17.711 . Note that for each curve, the values are arranged in ascending order so that, in general, 
a point on the horizontal axis does not correspond to the same target for all curves. 



30 AU (see their Appendix 2 for a discussion) and it 
is consistent with recen t results for stellar binaries be- 
tween ~ 5-500 AU fee. iKraus eFal. 2008. 201i:). Fur- 
thermore, in the event that >10Mj objects can form 
by the same mechanism as lower-mass gas giant plan- 
ets, Equa tion [T21 is a reasonable extrapolation from the 
results of iCumming et aTl (120081 ). who found a nearly- 
log-flat distribution for <10Mj gas giant planets at sep- 
arations a < 3 AU. 

7.6. Distribution-independent Approach 

In addition to assuming specific forms for the distribu- 
tion of companion properties, we repeated the analysis 
with them set to fixed values. This allowed us to obtain 
conservative upper limit estimates for the companion fre- 
quencies. For instance, the less massive a companion is, 
the more difficult it is to detect because it is fainter and 
hence the required contrasts are higher. Therefore, by 
setting all of our hypothetical substellar companions to 
have some mass iV/j, the subsequent result we obtain 
from the MC analysis is an upper limit on the frequency 
of all companions with masses M2 > Mj. 

Similarly, our ability to detect companions varied with 
angular separation (Figure [6]), which is related to the 
semimajor axis of the companion via the distance to the 
system and the orientation of the orbit. Now suppose we 
fix the semimajor axes of the hypothetical companions to 
a certain value a' and repeat the MC analysis for values 
over some interval a' S [oi , 02] . Then the maximum value 
of /„ obtained from these analyses is the most conserva- 
tive upper limit estimate for the frequency of companions 
with semimajor axes on that interval. 

We present the results of these distribution- 
independent calculations in Section 17. 9[ as well as 
the results obtained by assuming the specific distribu- 
tion forms described in Sections 17. 3H7.5I 

7.7. Previous Imaging Observations 

Ideally, when performing the calculations described 
above, we would like to combine the results of our aper- 



ture masking survey with those of other imaging sur- 
veys targeting wider angular scales. This would allow 
us to put tighter constraints on the companion frequen- 
cies across a larger range of separations. To this end, 
we identified 49 of our targets that have previously been 
observed as part of published direct imaging surveys and 
list these in Table [T^ For each of these targets, we quote 
the inner separation angle that was probed by the imag- 
ing observations as well as the corresponding sensitiv- 
ity of the observations. In most cases, these values were 
taken directly from the published survey limits, but when 
these were not provided explicitly we attempted to make 
conservative estimates. We also list each of the sensi- 
tivities in Table [T^l as an equivalent minimum detectable 
companion m ass, calculated by int erpolating the DUSTY 
isochrones of iChabrier et all (|2000 :) in the same manner 
outhned in Section 16.51 We incorp orat ed these limits 
into our analysis described in Section [721 by treating hy- 
pothetical companions as "detected" (i.e. by increasing 
Xj by 1) whenever they came within the detectability 
range of the previous imaging observations (i.e. if they 
had separation and masses above the values quoted in 
Table fT2|l. In the next sections, we present the results ob- 
tained from this combined approach (aperture masking 
-|- previous imaging) together with the results obtained 
using the aperture masking limits alone. 



7.8. Calculated pj Values 

The Pj values calculated separately for each of the 
three companion mass distributions that we considered 
(Equations O [H and [T0|) are plotted in ascending order 
in Figure [7] for the case of a companion orbital eccen- 
tricity distribution given by dAf/de oc 2e (Equation llip 
and semimajor axis distribution given by dAf/da oc 
(Equation [T2|) . In this figure, we immediately see the ad- 
vantage of combining our aperture masking results with 
the results from imaging surveys: the overall effect is 
roughly equivalent to an upwards shift of the pj values 
by ^ 10-30% for the majority of targets. 
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7.9. Calculated /„ Values 
7.9.1. Assuming dAf/da oc 

In Table [131 we present 95% confidence (i.e. a = 0.95 
in Equation [5]) upper limit estimates for the fre- 
quency of 20-80 Af J substellar companions in the sep- 
aration range 3-30 AU, with companion semimajor axes 
randomly drawn from the inverse power law distribution 
dJV/da oc a^^ (Equation [T2|). Also presented are cal- 
culations made separately for each permutation of the 
companion mass and eccentricity distributions described 
in Sections 17.31 and 1 7.41 respectively, as well as for differ- 
ent fixed companion masses of 20 Mj, 40 Mj, and 60 Mj 
(see Section FTG)) . 

All calculations reported in Table [13] are reasonably 
robust to the different assumptions made for the com- 
panion eccentricities, with the calculated upper limits 
only differing by <l-2% depending on whether all com- 
panion eccentricities are fixed to e = or e = 0.9, or if 
they are drawn randomly from a distribution of the form 
dAf/de oc 2e (Equation [TT]) . When the previous imaging 
observations are incorporated into the calculations and 
the ages and distances listed in Tables [5| and [31 are used, 
the upper limit estimates vary between 9-12%, depend- 
ing on which form is assumed for the distribution of com- 
panion masses, but irrespective of what is assumed about 
the orbital eccentricities. When the previous imaging ob- 
servations are not included in the analysis, the equiva- 
lent limits vary between 13-19%. For fixed companion 
masses of 20 Mj, 40 Mj and 60 Mj, the upper limit es- 
timates vary between 15-16%, 11-12% and 9-10%, re- 
spectively, when the imaging observations are included, 
and between 25-27%, 18-20% and 14-15%, respectively, 
when the imaging observations are not included. 

To investigate how sensitive our results are to the un- 
certainties in the distances and ages of our targets, we 
repeated the above calculations using the la upper lim- 
its for the ages and distances of each target provided in 
Tables [2] and [3] For instance, instead of using a distance 
of 28 pc and an age of 110 Myr for PW And, we used 
28 -f 7 = 35 pc and 110 4- 40 = 150 Myr, respectively. 
Assuming upper values for the ages and distances in this 
way results in a downward revision of our sensitivities to 
faint companions at smaller separations. Therefore, we 
had to re-calculate the survey detection limits presented 
in Table [TT] before repeating the analysis described in 
Sections 17.2117.71 Depending on which of the compan- 
ion mass distributions is used, the upper limit estimates 
obtained from this analysis vary between 11-15% when 
the imaging observations are included, and between 17- 
24% when the imaging observations are not included. 
For fixed companion masses of 20 Afj, 40 Afj and 60A'/j, 
when the imaging observations are included the calcu- 
lated upper limits vary between 20-23%, 14-15% and 
11-12%, respectively, and when the imaging observations 
are not included they vary between 38-40%, 24-25% and 
18-20%, respectively. 

We also investigated how sensitive our results are to 
the 9 targets of less certain membership identified in 
Section [3] (i.e. HD 89744, HD 92945, GJ466, EKDra, 
HIP 30030, TWA-21, TWA-6, TWA-14, TWA-23) by re- 
moving them and the 7 Her-Lyr targets from the anal- 
ysis. When all 16 of these targets are removed and we 
randomly sample the companion masses, the upper limit 



estimates vary between 12-15% when the imaging obser- 
vations are included and between 16-23% when the imag- 
ing observations are not included, depending on which of 
the three companion mass distributions from Section [7.31 
is used. For fixed companion masses of 20 A/j, 40 Mj and 
60 Afj, the upper limit estimates vary between 18-20%, 
14-15% and 12-13%, respectively, when the imaging ob- 
servations are included, and between 31-32%, 22-23% 
and 17-19%, respectively, when the imaging observations 
are not included. 

7.9.2. dM /da distribution-independent 

We repeated all of the calculations presented in the 
previous section for fixed companion semimajor axes be- 
tween 3-30 AU. As before, when the imaging observa- 
tions were included in the analysis, the upper limit esti- 
mates only change by ^1-2% over the entire 3-30 AU de- 
pending on which assumption is made for the companion 
eccentricities. However, when the aperture masking ob- 
servations are used on their own this variation increases 
to ^5-10% over the range 3-10 AU and becomes as high 
as -20% over the 10-30 AU range (Figure [S]). 

In Figure[S]we plot the calculated upper limit estimates 
obtained for fixed companion masses and randomly- 
sampled companion masses, while holding the semimajor 
axes fixed at successive values between 3-30 AU using a 
step size of 0.5 AU and randomly drawing the companion 
eccentricities from a distribution of the form dj\f/de oc 2e 
(Equation llip . On their own, the aperture masking re- 
sults place the tightest constraints over the — 3-lOAU 
semimajor axis range, with upper limit estimates of 20%, 
16% and 13% for fixed companion masses of 20Arj, 
40 M J and 60A/j, respectively. With the imaging ob- 
servations included in the analysis, these limits improve 
to 19%, 13% and 10%, respectively. At larger separa- 
tions between 10-30 AU, our upper limit estimates are 
12%, 9%, and 8%, respectively, for the same companion 
masses when we include the imaging observations, while 
the companion frequencies are poorly constrained by the 
aperture masking observations alone. 

Meanwhile, the right panel in Figure [H] shows the re- 
sults that were obtained when we sampled the companion 
masses from each of the three distributions given in Sec- 
tion [73] Over the 10-30 AU semimajor axis range, the 
aperture masking observations on their own constrain the 
frequency of 20-80 Afj companions to be less than 16%, 
15% or 13% at 95% confidence, depending on whether 
the mass power law (Equation [5]), mass ratio power law 
(Equation [5]) or mass ratio log- normal paramaterization 
(Equation [TP)) is assumed for the companions. These 
constraints improve to 13%, 12% and 10%, respectively, 
when the imaging observations are included in the anal- 
ysis. At wider separations between 10-30 AU, the equiv- 
alent values obtained when the aperture masking obser- 
vations are combined with the previous imaging observa- 
tions are 9%, 9% and 8%, respectively. 

Finally, Figures [TOl and [TT] have been included for com- 
pleteness. They are the same as Figure [H] except that 
they show respectively the results obtained when upper 
values for the target ages and distances are used as de- 
scribed in Section 17.9. 1[ and the results obtained when 
the 9 targets of less certain membership identified in Sec- 
tion [3] and the Her-Lyr targets are not included in the 
calculations. 
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Fig. 8. — 95% confidence fu estimates as a function of fixed companion semimajor axis obtained for the aperture masking limits alone 
with different assumptions about the companion eccentricities. The left panel shows values obtained for fixed companions masses of 20 Afj 
(black lines), 40Afj (red lines), and 60 Mj (green lines). The right panel shows the same values obtained for companions distributed 
according to the mass power law given by Equation [8] (bl ack lines), the mass ratio power law given by Equation [9] (red lines), and the 
log- normal mass ratio parameterization given by Equation 1101 (green lines). In both panels, fu values are shown for the following cases: 
fixed companion eccentricities of e = (solid lines); fixed companion eccen tric ities of e = 0.9 (dash-dot lines); and companion eccentricities 
randomly sampled from a distribution of the form dj\f /de oc2e (Equation [TT} (dashed lines). 



7.10. Implications for Formation Theories 

A well-known result from radial velocity surveys is 
the discovery of a "brown dwarf desert" at separations 
<3AU, where <0.5-l% of solar- like stars are found to 
possess a 13-75 Mj co mpanion (|Marcv fc Butleii [2OO0I : 
iGrether fc Lineweaveii[2006 ) com pared with '-^10% pos- 
sessing a 0.3-10 Mj companion (jCumming et al.l [20081) 
and ~13% possessing a > 0.1 Mcr, stellar companion 
(jDuquennov fc Mavojll99ll ). Meanwhile, at wider sep- 
arations, imaging surveys have started to place upper 
limits on the frequency of substellar companions: 



' iCarson et al.l (|2006f ) obtained a 95% confidence up- 
per limit of 12.1% on the frequency of 13-73 Afj 
companions between 25-100 AU. 



' iLafreniere et al.l (|2007D obtained a 95% confidence 
interval of 1.9l^;|% for the frequency of 13-75 Mj 
companions between 25-250 AU. 



• iMetchev fc Hillenbrandl (|2009[ ) obtained a 95% 
confidence interval of 3.2^2 7% for the frequency of 
13-75 Mj companions between 28-1590 AU This i s 
consistent with the results of iKraus et al.l ()2011[ ). 
who obtained a lower bound of 3.9^12% for the fre- 
quency of substellar companions over the range 5- 
5000 AU by combining the results of their aperture 
masking survey of Taurus- Auriga members with 
previous direct imaging results. 

• By jointly analyzing the results from three of the 
largest and deepest surveys for substellar com- 



panio n s to date (.Masciadri et al. 2005.; Biller et al 
[2007t ILafreniere et all l2007t ). iNielsen fc Close 
(|2010D obtained 95% confidence upper limits of 
<20% and <5% for the frequency of companions 
with masses between 10-15 Mj in the ranges 13- 
600 AU and 40-200 AU, respectively 



The aperture masking survey reported in this paper has 
allowed us to place similar constraints on the frequency of 
20-80 Mj companions over the 3-30 AU separation range 

(Sections [LSHISl)- 

These results are broadly in line with expectations 
from current models of substellar companion formation. 
Firstly, population synthesis models predict that core 
accretion only produces companions with masses up to 
-10 Mj (|Ida fc Linll200l . or else, if objects are formed 
with masses above 20 M j , then these are extremely rare 
(jMordasini et al.l 120091 ). Unsurprisingly, the observa- 
tional studies outlined above provide no evidence to the 
contrary, despite the aperture masking surveys in partic- 
ular probing the range of separations where core accre- 
tion is expected to be most efficient. 

Indeed, 20-80 Mj companions are much more likely to 
form by either gravoturbulent fragmentation during the 
initial collapse of the molecular cloud (eg. IBatd 120091 : 
lOffner et al.l 120091 ) or by the fragmentation of gravita- 
tional instabilities in the protostellar disk once the ini- 
tial free-fall collapse of the molecular cloud has ended 
(eg. IClarkd[2009l : IStamatellos fc WhitwOTth|[2009l) . For 
the gravoturbulent fragmentation scenario, the low fre- 
quencies of substellar companions deduced for separa- 
tions <200AU from observational studies is in quali- 
tative agreement with the hydrodynamical simulations 
of Bate (2009) who found that the separation of binary 
pairs consisting of a stellar primary and a very low-mass 
secondary increases strongly with decreasing mass ratio. 
Meanwhile, the disk fragmentation mechanism is not ex- 
pected to occur within ~40-70 AU of the primary, where 
radiative cooling time scale s are too long for the dis k 
to be Toomre unstable (eg. iRafikovl [2007t iBolevI 12009) . 
Alternatively, 20-80 Mj objects might form by gravita- 
tional disk instabilities at sepa rations beyond ~40-70 AU 
and t hen migrate inwards. IStamatellos fc WhitworthI 
(j2009| ) considered this for the case of a massive disk ex- 
tending between 40-400 AU around a 0.7 M0 star. How- 
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Fig. 9. — 95% confidence /„ estimates for companions with masses in the range 20-80 A/j as a function of fixed companion semimajor 
axes. As in Figure [7] the left panel shows the results for fixed companion masses and the right panel shows the results for companion 
masses drawn randomly from distributions. Colors and linestyles are the same as Figure [7] 
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Fig. 10. — The same as Figure[9] but using values for the target distances and ages that give conservative estimates for the companion 
detection sensitivities. Specifically, the upper la limits given in Tables [2] and |3] were used for the target ages and distances, respectively. 




Fig. 11. — The same as Figure |9] but excluding 16 of the ta rgets whose moving group membership is not well-established; namely, 
the 9 targets that did not meet the 'high probability' criterion of lTorres et al.l |200a ) in their dynamical convergence analysis (HD 89744, 
HD 92945, GJ 466, EK Dra, HIP 30030, TWA-21, TWA-6, TWA-14, TWA-23), and the 7 targets belonging to the proposed Her-Lyr moving 
group (see Section [St. 
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ever, they found that when low mass (<80Mj) com- 
panions did form at closer separations, they were sub- 
sequently scattered outwards by dynamical interactions 
with more massive companions in the same disk, lead- 
ing to a brown dwarf desert that extended out to ~100- 
200 AU. Again, the low occurence of 20-80 A/j compan- 
ions inferred from observational studies over this separa- 
tion range is consistent with such a scenario, though the 
constraints are not yet tight enough to make a definitive 
statement. 

8. CONCLUSION 

This paper has presented the results of an aperture 
masking survey of 67 young nearby stars for substel- 
lar companions. Our detection limits extend down to 

40 Af J for 30 of our targets, and of these, we are sen- 
sitive down to '^20A/j or less for a subset of 22. Al- 
though we did not uncover any substellar companions, 
we detected four stellar companions. One of these, a 
0.52 ± 0.09 Mq companion to HIP 14807, is a new discov- 
ery. We have also shown that the companion to HD 16760 
is on a low inclination orbit with a mass of O.28±O.O4M0, 
much higher than the minimum mass of AI2 sini ~13- 
14A/j inferred from radial velocity measurements. 

If we do not make any assumptions about the distri- 
bution of companion masses or semimajor axes, we cal- 
culate that the frequency of 20-80 Af j companions is less 
than ~ 19% in the range 3-10 AU and less than ^ 12% in 
the range 10-30 AU at 95% confidence. If, however, we 
assume that the semimajor axes of 20-80 A/j companions 
are distributed according to dM /da oc and that their 
masses are distributed according to a log-normal param- 



eterization of the secondary-to-primary mass ratio, this 
limit becomes ~9% over the 3-30 AU separation range. 
Similar values of ~10% and ~11% are obtained if we as- 
sume instead that the companion masses or secondary- 
to-primary mass ratios, respectively, are distributed ac- 
cording to power laws. These results are consistent with 
models that predict a low occurrence of substellar com- 
panions relative to stellar companions at these separa- 
tions, possibly hinting at the extension of the brown 
dwarf desert beyond ~3AU. 
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TABLE 1 
Mask Properties 



Observatory Telescope Instrument Mask Hole Diameter Transmission 



Diameter (m) (m) % 



Palomar 5.1 PHARO 9H 0.4 6 

18H 0.2 3 

Keck 10 NIRC2 9H 1.1 11 

18H 0.5 5 

VLT 8.2 CONICA 7H 1.2 15 



TABLE 2 
Moving Groups Surveyed 



Group 




Age^ 


Distance 


Reference 






(Myr) 


(pc) 




AB Dor 


31 


110 ±40 


34.1 ± 12.8 


Lu05, T08 


13 Pic 


11 


12 ±5 


34.5 ± 1.4 


ZOla, F06, T08 


Her-Lyr 


7 


200 ± 80 


14.6 ±4.1 


LS06 


Tuc-Hor 


2 


30 ± 10 


45.1 ± 0.6 


ZOlb, T08 


TWA 


16 


8±4 


55.0 ±2.7 


dR06, T08 



^ The number of targets observed for each group in this survey. 
^ Adopted age uncertainties are estimated from the spread in 
reported hterature values and their associated uncertainties. 
^ Quoted values and dispersions are the medians and standard 
deviations of those members observed in this survey (see Ta- 
ble [3} . 

^ ITorres et all ( 120081) (T08) provide a review of the 
memberships of four of the five groups targeted in this 
survey, while iLopez- Santiago et al.l 1 120061) (LS06) describe 
the proposed Her-Lyr association. Also cited here are 
Luliman et al.l ll2005ll fLuOSl. lde la Reza et all ll2006l) (dR06) , 
Zuckerman et al . |2b01a) ( ZOla) , IZuckerman et all j2001bD 
(ZOlb), and lFeigclson ct all | [20Ml (F06). 



TABLE 3 
Sample List 



03 



Name 


RA 


DEC 


Spectral 


Mass'' 


Distance 


Ref.'^ 


K 


H 


J 


Ref.d 


Membership 


Previous 




(J2000) 


Type^ 


(Mo) 


(pc) 




(mag) 


(mag) 


(mag) 




Ref.« 


Imaging * 


AB Dor Targets 


PW AnH 


on 1 s 90 8 

UU -LO ^KJ.O 




K2 


0.81 


28 + 78 


ZS04 


6.39 


6.51 


7.02 


2M 


7,9n4 TC\R 

ZiOU^5 -L UO 


M7,n4 T,n^ T,n7 
























MH09, HIO 


HTP S'iRQ 


nn 45 50 Q 


-i-U i^s 4(1 


F8 


1.12 


CO 1^ + 2 5 


HIP 


6.69 


6.72 


6.93 


2M 


ZSn4 T08 


HIP 5191 


01 06 26.1 


— 14 17 46 


Kl 


0.86 


47.3 ± 2.8 


HIP 


7.34 


7.43 


7.91 


2M 


ZS04, T08 


CIO 


HIP 6276 


01 20 32.2 


-11 28 03 


G9 


0.89 


34.4 ± 1.2 


HIP 


6.55 


6.65 


7.03 


2M 


ZS04, T08 


MH09 


HIP 12635 


02 42 20.9 


+38 37 22 


K2 


0.79 


50.4 ± 6.7 


HIP 


7.76 


7.90 


8.38 


2M 


ZS04, T08 




HD 16760 


02 42 21.3 


+38 37 08 


G5 


0.91 


45.5 ± 4.9 


HIP 


7.11 


7.15 


7.47 


2M 


ZS04, T08 




HIP 13027 


02 47 27.4 


+19 22 19 


Gl 


1.02 


33.6 ± 0.9 


HIP 


6.05 


6.10 


6.37 


2M 


ZS04, T08 




HD 19668 


03 09 42.3 


-09 34 47 


GO 


0.90 


37.4 ± 1.6 


HIP 


6.70 


6.79 


7.16 


2M 


LS06, T08 


MH09 


HIP 14807 


03 11 12.3 


+22 25 24 


K6 


0.76 


52.5 ± 8.6 


HIP 


7.96 


8.10 


8.67 


2M 


ZS04, T08 




HIP 14809 


03 11 13.8 


+22 24 58 


G5 


1.04 


53.7 ± 3.3 


HIP 


6.97 


7.07 


7.27 


2M 


ZS04, T08 




HIP 16563A 


03 33 13.5 


+46 15 27 


G5 


0.88 


34.4 ± 1.2 


HIP 


6.62 


6.70 


7.03 


2M 


ZS04, T08 


B07 


HIP 16563B 


03 33 14.0 


+46 15 19 


MO 


0.55 


34.4 ± 1.2 


HIP 


8.07 


8.21 


8.83 


2M 


ZS04, T08 




HIP 17695 


03 47 23.2 


-01 58 18 


M3 


0.46 


16.1 ± 0.7 


HIP 


6.93 


7.17 


7.80 


2M 


ZS04, T08 


L07 


HIP 18859 


04 02 36.7 


-00 16 06 


F6 


1.16 


18.8 ± 0.1 


HIP 


4.18 


4.34 


4.71 


2M 


ZS04, T08 


L07, HIO 


HIP 19183 


04 06 41.5 


+01 41 03 


F5 


1.17 


55.2 ± 2.8 


HIP 


6.58 


6.70 


6.89 


2M 


ZS04, T08 


BD+20 1790 


07 23 44.0 


+20 25 06 


K5 


0.76 


32 ± 88 


LS06 


6.88 


7.03 


7.64 


2M 


LS06, T08 


L05, MH09, HIO 


HD 89744 ' 


10 22 10.6 


+41 13 46 


F7 


1.52 


39.4 ± 0.5 


HIP 


4.45 


4.53 


4.86 


2M 


LS06 




HIP 51317 


10 28 55.6 


+00 50 28 


M2 


0.43 


7.1 ± 0.1 


HIP 


5.31 


5.61 


6.18 


2M 


LS06, T08 


M05, L07 


HD 92945'' 


10 43 28.3 


-29 03 51 


Kl 


0.85 


21.4 ± 0.3 


HIP 


5.66 


5.77 


6.18 


2M 


LS06 


B07, L07 


GJ466'' 


12 25 58.6 


+08 03 44 


MO 


0.73 


37.4 ± 3.2 


HIP 


7.33 


7.45 


8.12 


2M 


LS06 


MZ04 


HD 113449 


13 03 49.8 


-05 09 41 


Kl 


0.84 


21.7 ± 0.4 


HIP 


5.72 


5.89 


6.27 


2M 


ZS04, T08 


L07, HIO 


EK Dral" 


14 39 00.2 


+64 17 30 


G1.5 


1.06 


34.1 ± 0.4 


HIP 


5.91 


6.01 


6.32 


2M 


LS06 


MZ04, B07, L07, 
MH09 


HIP 81084 


16 33 41.7 


-09 ,33 10 


MO 


0.58 


30.7 ± 2.3 


HIP 


7.55 


7.78 


8.38 


2M 


ZS04, T08 


L07 


HIP 82688 


16 54 08.2 


-04 20 24 


GO 


1.12 


46.7 ± 2.0 


HIP 


6..36 


6.48 


6.70 


2M 


ZS04, T08 


MH09 


HD 160934 


17 38 39.7 


+61 14 16 


K7 


0.70 


33.1 ± 2.2 


HIP 


7.22 


7.37 


7.98 


2M 


ZS04, T08 


L05, H07, L07, 
MZ04 


HIP 106231 


21 31 01.6 


+23 20 09 


K5 


0.75 


24.8 ± 0.7 


HIP 


6.38 


6.52 


7.08 


2M 


ZS04, T08 


LOS, L07, MZ04 


HIP 110526 


22 23 29.1 


+32 27 34 


M3 


0.48 


15.5 ± 1.6 


HIP 


6.05 


6.28 


6.90 


2M 


ZS04, T08 




HIP 113579 


23 00 19.2 


-26 09 12 


G5 


0.99 


30.8 ± 0.7 


HIP 


5.94 


6.04 


6.29 


2M 


ZS04, T08 


MH09, CIO 


HIP 114066 


23 06 04.6 


+63 55 35 


Ml 


0.60 


24.5 ± 1.0 


HIP 


6.98 


7.17 


7.82 


2M 


ZS04, T08 


L07 


HIP 115162 


23 19 39.5 


+42 15 10 


G4 


0.94 


50.2 ± 2.9 


HIP 


7.22 


7.28 


7.61 


2M 


ZS04, T08 




HIP 118008 


23 56 10.5 


-39 03 07 


K2 


0.81 


22.0 ± 0.4 


HIP 


5.91 


6.01 


6.51 


2M 


ZS04, T08 


B07, CIO 



13 Pic Targets 



HR9 


00 


06 


50, 


.1 


-23 


06 


27 


F3 


1, 


.40 


39, 


,4 


± 


0, 


,6 


HIP 


5. 


,24 


5. 


,33 


5.45 


2M 


ZS04, 


T08 


K03 


HIP 10680 


02 


17 


25, 


.2 


+28 


44 


43 


F5 


1, 


,15 


34, 


,5 


± 


0, 


,6 


HIP 


5. 


,79 


5, 


,84 


6.05 


2M 


ZS04, 


T08 




HIP11437B 


02 


27 


28, 


,1 


+30 


58 


41 


M2 


0, 


,35 


40, 


,0 


± 


3, 


,6 


HIP 


7, 


,92 


8, 


,14 


8.82 


2M 


ZS04, 


T08 




HIP11437A 


02 


27 


29, 


,2 


+30 


58 


25 


K6 


0, 


,63 


40, 


,0 


± 


3, 


,6 


HIP 


7, 


,08 


7, 


,24 


7.87 


2M 


ZS04, 


T08 




HIP 12545 


02 


41 


25, 


,8 


+05 


59 


19 


K6 


0, 


,67 


42, 


,0 


± 


2, 


,7 


HIP 


7, 


,07 


7, 


,23 


7.9 


2M 


ZS04, 


T08 


B07 


51 Eri 


04 


37 


36, 


,1 


-02 


28 


25 


FO 


1, 


,41 


29, 


,4 


± 


0, 


,3 


HIP 


4, 


,54 


4, 


,77 


4.74 


2M 


ZS04, 


T08 


HIO 


HIP 25486 


05 


27 


04, 


,8 


-11 


54 


04 


F7 


1, 


,25 


27, 


,0 


± 


0, 


,4 


HIP 


4. 


,93 


5, 


,09 


5.27 


2M 


ZS04, 


T08 


L05, K07, MH09 


GJ 803 


20 


45 


09, 


.5 


-31 


20 


27 


Ml 


0, 


.44 


9, 


,9 


± 


0, 


,1 


HIP 


4. 


,53 


4, 


,83 


5.44 


2M 


ZS04, 


T08 


K03, MZ04, M05, 




















































B07, L07 


BD-17 6128 


20 


56 


02, 


,7 


-17 


10 


54 


K6 


0, 


,73 


45, 


,7 


± 


1, 


,6 


HIP 


7. 


12 


7, 


,25 


7.92 


K04 


ZS04. 


T08 


M05 


HIP 112312A 


22 


44 


57, 


,9 


-33 


15 


02 


M4 


0, 


,31 


23, 


,3 


± 


2, 


,0 


HIP 


6. 


,93 


7, 


,15 


7.79 


2M 


ZS04, 


T08 


B07 


HIP 112312B 


22 


45 


00, 


,0 


-33 


15 


26 


M4.5 


0, 


,17 


23, 


,3 


± 


2, 


,0 


HIP 


7. 


79 


8, 


,06 


8.68 


2M 


ZS04, 


T08 





Her-Lyr Target! 



TABLE 3 
Sample List 



HD166 


00 06 36.8 


29 01 17.4 


KO 


0.93 


13.7 ± 0.1 


HIP 


4.31 


4.63 


4.73 


2M 


LS06 


L07, HIO 


HD 10008 


01 37 35.5 


-06 45 37.5 


G5 


0.89 


24.0 ± 0.4 


HIP 


5.75 


5.90 


6.23 


2M 


LS06 


L07 


HD 233153 


05 41 30.7 


+53 29 23 


MO. 5 


0.58 


12.4 ± 0.3 


HIP 


5.76 


5.96 


6.59 


2M 


LS06 


C05 


xllr^ o ( zoo 


(T? Qn OQ n 


1 no 1 1 ni 


T<"7 
IV / 


U.Ol 


14.0 it U.o 






a no 


Ki.il 


ZiVl 


L-oUO 


iViUO, Li(J 1 


HD 70573 


08 22 50.0 


01 51 33.6 


G6 


0.89 


46 ± lis 


LS06 


7.19 


7.28 


7.56 


2M 


LS06 


L05, MH09 


HIP 53020 


10 50 52.1 


+06 48 29 


M4 


0.25 


6.8 ± 0.2 


HIP 


6.37 


6.71 


7.32 


2M 


LS06 


L07 


HN Peg 


21 44 31.3 


+14 46 19 


GO 


1.06 


17.9 ± 0.1 


HIP 


4.56 


4.6 


4.79 


2M 


LS06 


MZ04, L07 


Tuc-Hor Targets 


























HIP 9141 


01 57 48.9 


-21 54 05 


G4 


0.97 


40.9 ± 1.1 


HIP 


6.47 


6.56 


6.86 


2M 


ZS04, T08 


B07, MH09 


HIP 30030 ^ 


06 19 08.1 


-03 26 20 


GO 


1.03 


49.2 ± 2.0 


HIP 


6.55 


6.59 


6.85 


2M 


ZS04 


B07, MH09 


TWA Targets 


























TWA-21 


10 13 14.8 


-52 30 54 


K3/4 


0.63 


48 ± 4 


MM05 


7.19 


7.35 


7.87 


2M 


ZS04 




TWA-6 


10 18 28.8 


-31 50 02 


K7 


0.43 


55 ± 5 


MM05 


8.04 


8.18 


8.87 


2M 


ZS04 


W99, MZ04, L05, 


TWA-7 


10 42 30.3 


-33 40 17 


M2 


0.35 


29 ± 2 


MM05 


6.9 


7.13 


7.79 


2M 


ZS04, T08 


M05 
W99, MZ04, L05 


TW Hya 


11 01 51.9 


-34 42 17 


K6 


0.64 


53.7 ± 6.2 


HIP 


7.30 


7.56 


8.22 


2M 


ZS04, T08 


W99, MZ04, L05 


TWA-3 


11 10 28.0 


-37 31 53 


M4 


0.37 


36 ± 4 


MM05 


7.28 


7.60 




WOO 


ZS04, T08 


W99, CIO 


TWA- 14 h 


11 13 26.5 


-45 23 43 


MO 


0.57 


86 ± 8 


MM05 


8.50 


8.73 


9.42 


2M 


ZS04 


B07, MZ04, CIO 


TWA-13B 


11 21 17.2 


-34 46 45 


Ml 


0.63 


57 ± 10 


MM05 


7.49 


7.73 


8.43 


2M 


ZS04, T08 




TWA-13A 


11 21 17.5 


-34 46 50 


Ml 


0.61 


57 ± 10 


MM05 


7.46 


7.68 


8.43 


2M 


ZS04, T08 




TWA-8B 


11 32 41.4 


-26 52 08 


M5 


0.14 


42 ± 5 


MM05 


9.01 


9.36 




WOO 


ZS04, T08 


W99, M05, L05 


TWA-8A 


11 32 41.5 


-26 51 55 


M3 


0.40 


41 ± 4 


MM05 


7.44 


7.72 




WOO 


ZS04, T08 


W99, M05 


TWA-9 


11 48 24.2 


-37 28 49 


K5 


0.38 


46.8 ± 5.4 


HIP 


7.85 


8.03 


8.68 


2M 


ZS04, T08 


W99, M05 


TWA-23'' 


12 07 27.4 


-32 47 00 


Ml 


0.58 


61 ± 5 


MM05 


7.75 


8.03 


8.62 


2M 


ZS04 


CIO 


TWA-25 


12 15 30.8 


-39 48 42 


Ml 


0.68 


55 ± 4 


MM05 


7.31 


7.50 


8.17 


2M 


ZS04, T08 


B07, CIO 


TWA-10 


12 35 04.3 


-41 36 39 


M2 


0.39 


57 ± 9 


MM05 


8.19 


8.48 


9.12 


2M 


ZS04, T08 


W99, MZ04, L05 


TWA-llB 


12 36 00.6 


-39 52 16 


M2 


0.52 


72.8 ± 1.7 


HIP 


8.35 


8.53 


9.15 


2M 


ZS04, T08 


W99 


TWA-llA 


12 36 01.0 


-39 52 10 


AO 


2.31 


72.8 ± 1.7 


HIP 


5.77 


5.79 


5.78 


2M 


ZS04, T08 


W99, CIO 



TABLE 3 
Sample List 



00 



Spectral types arc taken from the lists contained in lZuckerman &: Sond i2004D . IL6pez- Santiago et all i2Q06f) . and I Torres et all i200SD . 

Masses are estimated by interpolation of the NextGen isochrones of lBaraffe et ahl j 19981) , assuming the ages listed in Table[2] The only exeeptions are HD 89744 and TWA-llA whieh fall 
outside the range of the NextGen isochrones. For these two targets, we estimated their masses using the Y"^ isochrones of iYi et all (]200in . Allowing for the uncertainties in the distances and 
ages, the calculated masses are typically uncertain at the <5— 10% level for the AB Dor targets, ^30% for the /3 Pic targets, ^5% for the Her-Lyr and Tuc-Hor targets, and <40% for the TWA 
targets. These values, however, do not account for uncertainties in the isochrones themselves, nor for uncorrected blending. 

^ Distances have been obtained from HIPPARCOS (HIP) parallaxes where available jvan Leeuwen||2007D . In the case of BD-17 6128, the HIPPARCOS distance to HD 199143 is quoted, 
since [van den Ancker et all (]2000l ) have shown that they form a physical pair. Otherwise distances have been derived photometrically, as listed bv iZuckerman fc Son^ <]2004|) (ZS04), 
iLopez-Santiago et al.| d2006f) (LS06), and iMamajekl ||2005D (MM05). Note that here the MM05 distances have been increased upwards by 7% following the revision of iMamajek & Meyer! 
(f2007l) . 

2MASS (2M) photometry is used where available jCutri et al.|[2003r . A blend correction is applied in cases where the target is known to be binary but is unresolved in the 2MASS catalog. 
For other targets, quoted magnitudes were obtained from lKaisler et al.| j2004D (K04) and IWeintraub et all (120001) (WOO). 

® IZuckerman & Sond j2004D (ZS04). [TTopez- Santiago et al.| j2006D (LS06), and ITorres et al.| (120081) . A number of targets were originally on the lists of ZS04 or proposed by LS06, but not 

confirmed by T08. 
f 

Previous high-resolution imaging surveys that have observed the target. Most used coronagraphic or differential imaging techniques to search for low-mass companions at wider separations that 
are complementary to ours. IBiller et all j2007D fB07) . fCarson et al.| ( f2005l) fC05). fChauvin et al.] ( i20Tol) fClO). [Hormuth et all ( f2007l) (H07),lHeinzc et al.' ('2010b||a|) (HIO) . [Kaisler et al.| ( i2003l ) 
('K03), fKasper et al.| j2007D f K07) . [Lowrance et ahl | |2005D (LOS). [Lafreniere et all I I2007D ('L07). [Masciadri et all (120051) fM05). [Metchev & HTTlcnbrand L200a) fMH09). [McCarthv & ZuckermanI 
( f200l) fMZ04). |mbb et al.| | fT99g|) (W99). 

^ Uncertainty of 25% assumed for photometrically-determined distances. 

^ Proposed members that were not confirmed by the dynamical convergence method of T08. The analysis described in Section 17. 2| was repeated separately both with and without these 9 
targets included. ^ 
^ HD 89744 is orbited by a substellar companion (Af2sinz ~7— 8Afj) with a 256 day period jKorzennik et al.l|2000l : iButler et al.l|2006D . Isochrone fits suggest that the system is ~2 Gyr old ^ 
(]Ng &: Bertellll99Sl : [Gonzalez et al.l|2001l) . and hence does not belong to the AB Dor moving group, as suggested by LSQ6. This older age is supported by evidence from gyrochronology jBarne j ^ 
120071) . Furthermore, T08 did not find strong evidence for HD 89744 being a member of AB Dor on dynamical grounds. However, we have chosen to retain HD 89744 in our sample for two ^ 
reasons: (1) with a scmimajor axis of 0.88 AU and a distance of 39.4 pc, the companion was potentially within our range of detectability; however, (2) in the end we were only sensitive to p 
companions with masses >200Mj around HD 89744 ( Table [TTl l. and as such, its inclusion does not affect the calculations presented in Scction[7j as these only included targets with sensitivity 
limits extending below 80 Mj . 
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TABLE 4 

Filter H UsEn 



Instruincrit 


Filter 


A() (/im) 


AA ((Um) 


PHARO 


CH4s 


1.57 


0.10 




H 


1.64 


0.30 




Ks 


2.15 


0.31 


NIRC2 


Jcont 


1.21 


0.02 




Hcont 


1.58 


0.02 




CH4s 


1.59 


0.13 




Kp 


2.12 


0.35 




Kcoiit 


2.27 


0.03 




CO 


2.29 


0.03 


CONICA 


L' 


3.80 


0.62 



TABLE 5 
Observations 



Target Instrument 


Filter 


Mask 


Tint (min) 


Date'' 




MJD 


AB Dor Targets 


PW And 


PHARO 


CH4s 


18H 


3.9 


2007 Jun 


01 


54252, 


,5 


HIP 3589 


PHARO 


CH4s 


9H 


6.5 


2007 Nov 


27 


54431, 


,1 


HIP 5191 


NIRC2 


Hcont 


18H 


4.0 


2007 Jun 


06 


54257, 


,6 


HIP 6276 


NIRC2 


Kcont 


18H 


5.3 


2007 Nov 


23 


54427, 


,2 


HIP 12635 


NIRC2 


Kcont 


18H 


800 


2008 Dec 


23 


54823, 


,2 


HD 16760*' 


NIRC2 


Kcont 


18H 


2.7 


2008 Dec 


23 


54823, 


,2 




NIRC2 


Kcont 


18H 


13 


2009 Aug 


06 


55049, 


.6 




NIRC2 


Jcont 


18H 


2.7 


2009 Nov 


20 


55155, 


,3 




NIRC2 


Hcont 


18H 


2.7 


2009 Nov 


20 


55155, 


,3 




NIRC2 


Kcont 


18H 


2.7 


2009 Nov 


20 


55155, 


.3 




NIRC2 


CO 


18H 


2.7 


2009 Nov 


21 


55156, 


.2 


HIP 13027 


PHARO 


Ks 


9H 


19 


2007 Nov 


27 


54431, 


,2 


HD 19668 


PHARO 


Ks 


9H 


5.8 


2007 Nov 


29 


54433, 


,3 


HIP 14807*= 


PHARO 


Ks 


9H 


5.8 


2007 Nov 


29 


54433, 


,1 




NIRC2 


CO 


18H 


5.0 


2009 Nov 


21 


55156, 


,2 


HIP 14809 


NIRC2 


CO 


18H 


2.7 


2009 Nov 


21 


55156, 


,2 


HIP 16563A 


PHARO 


Ks 


9H 


5.8 


2007 Nov 


27 


54431, 


,3 


HIP 16563B 


NIRC2 


CO 


18H 


5.3 


2009 Nov 


21 


55156, 


,2 


HIP 17695 


PHARO 


Ks 


9H 


12 


2007 Nov 


29 


54433, 


,3 


HIP 18859 


PHARO 


Ks 


9H 


8.7 


2007 Nov 


27 


54431, 


,2 


HIP 19183 


PHARO 


Ks 


9H 


8.7 


2007 Nov 


27 


.54431, 


,3 


BD+20 1790 


PHARO 


Ks 


9H 


5.8 


2007 Nov 


27 


54431, 


,5 


HD 89744 


PHARO 


CH4s 


18H 


8.6 


2007 Apr 


05 


54195, 


.3 


HIP 51317 


PHARO 


Ks 


9H 


13 


2007 Apr 


06 


54196, 


.3 


HD 92945 


NIRC2 


Kcont 


18H 


1.3 


2007 Jun 


06 


54257, 


.2 


GJ466 


PHARO 


Ks 


9H 


2.2 


2008 Jun 


19 


54636, 


.2 


HD 113449*' 


PHARO 


CH4s 


18H 


13 


2007 Apr 


06 


54196, 


.3 




PHARO 


Ks 


9H 


7.5 


2007 Apr 


07 


54197, 


.4 




PHARO 


CH4s 


18H 


3.9 


2007 Jun 


01 


54252, 


,1 




NIRC2 


Hcont 


18H 


1.3 


2008 Jun 


17 


54634, 


,2 




NIRC2 


CH4s 


9H 


1.7 


2008 Dec 


21 


54821, 


.7 




NIRC2 


Hcont 


18H 


2.7 


2010 Apr 


25 


55311, 


.4 



20 



Evans et al 
TABLE 5 — Continued 



Target 


Instrument 


Filter 


Mask 


Ti„t (min) 


Date'' 


MJD 


EK Dra 


PHARO 


Ks 


9H 


4.3 


2008 Jun 20 


54637.2 


HIP 81084 


PHARO 


CH4s 


18H 


4.3 


2007 May 30 


54250.2 


HIP 82688 


PHARO 


CH4s 


18H 


4.3 


2007 Jun 02 


54253.4 


HD 160934'' 


PHARO 


H 


9H 


9.7 


2008 .Jun 23 


54640.3 




PHARO 


Ks 


9H 


8.6 


2008 .Jun 23 


54640.3 




NIRC2 


Kcont 


Clear 


0.2 


2010 Apr 26 


55312.6 




NIRC2 


Jcont 


18H 


2.7 


2011 Apr 23 


55674.6 




NIRC2 


Hcont 


18H 


2.7 


2011 Apr 23 


55674.6 


HIP 106231 


PHARO 


CH4s 


18H 


4.8 


2007 May 31 


54251.4 


HIP 110526 


PHARO 


CH4s 


9H 


1.9 


2007 May 31 


54251.4 


HIP 113579 


NIRC2 


Hcont 


18H 


2.7 


2007 Jun 05 


54256.6 


HIP 114066 


PHARO 


CH4s 


18H 


3.9 


2007 Jun 01 


54252.4 


HIP 115162 


NIRC2 


Kcont 


18H 


5.3 


2009 Nov 21 


55156.2 


HIP 118008 


NIRC2 


Kcont 


18H 


5.3 


2007 Nov 23 


54427.2 


/3 Pic Targets 














HR9 


NIRC2 


Hcont 


18H 


2.7 


2007 Jun 06 


54257.6 


HIP 10680 


NIRC2 


Kcont 


18H 


5.3 


2007 Nov 24 


54428.4 


HIP 11437A 


NIRC2 


Kcont 


18H 


5.3 


2007 Nov 24 


54428.4 


HIP 11437B 


NIRC2 


Kcont 


18H 


5.3 


2007 Nov 24 


54428.4 


HIP 12545 


NIRC2 


Kcont 


18H 


5.3 


2007 Nov 24 


54428.4 


51 Eri 


NIRC2 


Kcont 


18H 


5.3 


2008 Dec 21 


54821.4 


HIP 25486 


PHARO 


Ks 


9H 


12 


2007 Nov 27 


54431.4 


GJ 803 


NIRC2 


Hcont 


18H 


4.0 


2007 Jun 05 


54256.5 


BD-17 6128 


PHARO 


CH4s 


9H 


4.3 


2007 May 30 


54250.5 


HIP 112312A 


NIRC2 


Kp 


9H 


10 


2008 Jun 17 


54634.6 


HIP 112312B 


NIRC2 


Hcont 


18H 


2.7 


2007 Jun 05 


54256.6 


rler-Lyr targets 














HD166 


PHARO 


CH4s 


18H 


5.8 


2007 May 31 


54251.5 


HD 10008 


NIRC2 


Kcont 


18H 


5.3 


2007 Nov 23 


54427.3 


HD 233153 


PHARO 


Ks 


9H 


8.7 


2007 Nov 27 


54431.5 


HIP o72o8 




Kcont 


18H 


5.3 


2007 JNov 24 


54428.0 


HD 70573 


PHARO 


Ks 


9H 


13 


2007 Nov 27 


54431.5 


HIP 53020 


NIRC2 


Kcont 


18H 


2.3 


2007 Jun 06 


54257.3 


HN Peg 


PHARO 


CH4s 


18H 


3.9 


2007 Jun 01 


54252.5 


Tuc-Hor Target! 


3 












HIP 9141 


PHARO 


Ks 


9H 


9.7 


2007 Nov 29 


54433.2 


HIP 30030 


NIRC2 


Kcont 


18H 


5.3 


2007 Nov 24 


54428.6 


TWA Targets 














TWA-21 


CONICA 


L' 


7H 


20 


2009 Mar 06 


54896.1 


TWA-6 


CONICA 


L' 


7H 


20 


2009 Mar 07 


54897.1 


TWA-7 


CONICA 


L' 


7H 


40 


2009 Mar 06 


54896.1 


TW Hya 


CONICA 


V 


7H 


60 


2009 Mar 07 


54897.2 


TWA-3 


NIRC2 


Kp 


9H 


11 


2008 Dec 22 


54822.7 


TWA- 14 


CONICA 


L' 


7H 


20 


2009 Mar 07 


54897.1 


TWA-13B 


NIRC2 


Kp 


9H 


5.3 


2008 Dec 21 


54821.6 


TWA- 13 A 


CONICA 


L' 


7H 


40 


2009 Mar 06 


54896.2 


TWA-8B 


NIRC2 


Kp 


9H 


5.3 


2008 Dec 23 


54823.7 


TWA-8A 


CONICA 


L' 


7H 


20 


2009 Mar 05 


54895.3 


TWA-9 


CONICA 


L' 


7H 


20 


2009 Mar 05 


54895.3 


TWA-23 


CONICA 


L' 


7H 


25 


2009 Mar 05 


54895.3 


TWA-25 


CONICA 


L' 


7H 


20 


2009 Mar 06 


54896.2 


TWA-10 


CONICA 


L' 


7H 


20 


2009 Mar 06 


54896.3 


TWA-llB 


CONICA 


L' 


7H 


20 


2009 Mar 07 


54897.3 


TWA-llA 


CONICA 


L' 


7H 


50 


2009 Mar 07 


54897.4 



^ For those targets that were observed on multiple epochs, we only report the observation that gave the deepest limits. The only exceptions are 
for those targets with detected companions, for which all epochs are reported. 
^ Companion detected in the current survey. 

Observation made in imaging mode without aperture mask. 



TABLE 6 
Detected Companions 



Primary 


Instrument 


Date 


MJD 


Filter 


Ao (^m) 


p (mas) 


Q (dcg) 


C 




TTTTll OnT 


T>1U A D r\ 


on/TV "XT „ or» 

z(JU7 INov z9 


o44oo.l 


Ks 


2.15 


63.22 


246.86 


10.15 






NIRC2 


2009 Nov 21 


55156.2 


CO 


2.29 


28.74 ±0.19 


89.74 ±0.29 


3.00 ±0.06 


0.52 ±0.09 


TTT~V1 fTrTi 

HDIdydU 




zOOo Dec 2o 


C A OOO O 

54823.2 


Kcont 


2.27 


o/^ 11 1 o cm 

2o.ll ± 2.59 


4d.2U ± 1.2d 


1 1 /I O 1 O 1 CT 

13.48 ± 3.15 


c\ lo 1 nil 

U.32 ± U.ll 






zUUy Aug D 


OoU4y.u 


Kcont 


O OT 


o^; 7Q _L n on 
2o. ( o iL u.yu 


zU4.a4 it U.4a 


1 Q 11 _L 1 nn 
iv5.il it l.UU 


n QO _L n no 

u.v>2 it u.uy 




NIRC2 


9nnQ Nn^r 90 


OO \-0\J . o 


J cont 


1.21 


OQ 1 "3 _|_ 1 QQ 


98*1 '^0 4- *i9 
zou.ou ^ o.uz 


Ql 04 _|_ t; qo 


n 94 -L 08 




NIRC2 


9nnQ Nn^r 90 


'J'-J i- oo .o 


r— r r^r\Y\\' 


1.58 


ofi nfi _|_ 1 jv. 


98*1 87+1 Q4 
zou . o / ^ ± . y± 


on CTQ _|_ 1 QQ 


fl 97 -U n 08 
u. z I ^ u.uo 






onnn IVT^,,^ on 
zUUy INov zU 


OOlOO.vS 


Kcont 


2.27 


39.07 


28D.D/ 


26.58 






NIRC2 


2009 Nov 21 


55156.2 


CO 


2.29 


25.37 ±4.51 


290.25 ±1.95 


15.18 ±6.44 


0.29 ±0.15 


TUTA 1 1 Q /I /I n 


Jr^rlAJrCU 


onnT A t-^v. R 
zUU/ Apr D 


r: -1 1 nt; q 


Url4s 


1.0 / 


or /?o _L n 1 
O0.D2 ± U.Ol 


ooc; in _L c\ A A 
22o. iy lb U.44 


A _L not; 
4.2 r ± U.2o 


n c; 1 _L n no 




"DTT A "D O 

Jr^rlAxxU 


onn'/ A r^t. v 
zUU/ Apr i 


o4iy / .4 


iVS 


lie; 


/in /If; 
4U.4u 


00*3 Qn 


U.DO 






"OTT A ID r\ 


zOUy Jun 1 


/I o cr o 1 b 

54252.1 


CH4s 


1.57 


28.46 


179.92 


13.73 








onno 1 '7 

zOUo Jun 17 


54do4.2 


Hcont 


1.58 


or" ^ o 1 n 1 o 

do. 68 ± 


222.91 ± U.21 


/I c CT 1 n n cr 

4.65 ± 0.U5 


n crn i n no 

U.5U ± U.02 




NIRC2 


2008 Dec 21 


54821.7 


CH4s 


1.59 


27.87 ±0.13 


250.14 ± 0.18 


4.62 ±0.03 


0.51 ±0.03 




NIRC2 


2010 Apr 25 


55311.4 


Hcont 


1.58 


35.81 ±0.17 


202.38 ± 0.23 


4.58 ±0.06 


0.51 ±0.02 


HD160934 


PHARO 


2008 Jun 23 


54640.3 


H 


1.64 


169.24 ±0.13 


273.35 ± 0.05 


2.22 ±0.01 


0.54 ±0.03 




PHARO 


2008 Jun 23 


54640.3 


Ks 


2.15 


169.79 ±0.25 


273.29 ± 0.09 


2.11 ±0.02 


0.54 ±0.04 




NIRC2 


2010 Apr 26 


55312.6^ 


Kcont 


2.27 


68.8 ±0.7 


290.0 ±0.6 


2.1 ±0.2 


0.54 ± 0.04 




NIRC2 


2011 Apr 23 


55674.6 


J cont 


1.21 


19.96 ±0.05 


18.44 ±0.12 


2.21 ±0.01 


0.54 ±0.03 




NIRC2 


2011 Apr 23 


55674.6 


Hcont 


1.58 


20.00 ± 0.03 


18.42 ± 0.09 


2.18 ±0.01 


0.54 ±0.03 



^ The quoted masses were calculated by interpolation of the NextGen models of lBaraffe et all (]1998D . using the primary magnitude in the approriate band as listed in Table [s] and the best-fit 
contrast C. Given that all of our companions were detected around AB Dor targets, we assumed a system age of 110 ± 40 Myr (Table [2j. Quoted uncertainties were calculated the same 
way as described in the footnotes of Table [s] The uncertainty associated with each mass value is the quadrature sum of three components: (1) AA/c* — \M{C ~ crc) — M{C + crc)|/2, 
where C is the fitted contrast with associated uncertainty ac', (2) AMy — \M{Y — (Ty) — M{Y + (7i-)|/2, where Y is the age with associated uncertainty cty as listed in Table [2] and (3) 
AAIfi — \M{d — ad) — M(d + (T(i)|/2, where d is the distance with associated uncertainty (T^ as listed in Column 6 of Table [3] 
^ Due to a contrast/separation degeneracy at close separations, these fits were unreliable and repeated with fixed contrasts (see Table [Tb. 

^ For this data set, the separation was well-determined but the contrast was poorly constrained due to degeneracy. The fit was repeated using the separation determined from the previous 
night (p ^ 35.62 ± 0.51mas) as a prior, giving a revised solution of C ^ 3.18 ± 1.59, p ^ 35.90 ± 0.47 mas and 6 ^ 223.09 ± 1.90 deg. The latter translates to a mass of M2 ^ 0.55 ± O.OSM© 
according to the NextGen isochrones of lBaraffe et all l |1998D . 

^ These data were taken with imaging, rather than aperture-masking, with the astrometry derived from the image auto-correlations. 
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TABLE 7 

Revised Fits For Degenerate Solutions 













Prior 




Refitted Solution 


Primary 


Instrument 


Date 


MJD 


Filter 


C 


C 


p (mas) 


e (dcg) 


HIP14807 


PHARO 


2007 Nov 29 


54433.1 


Ks 


3.00 ±0.06 


3.00 ±0.06 


45.43 ± 1.15 


248.08 ± 2.30 


HD16760 


NIRC2 


2009 Nov 20 


55155.3 


Kcont 


13.19 ±0.94 


14.05 ±0.92 


25.80 ± 1.03 


287.99 ± 1.89 


HD113449 


PHARO 


2007 Jun 01 


54252.1 


CH4s 


4.62 ±0.02 


4.62 ±0.02 


21.97 ± 0.73 


179.96 ±2.84 



TABLE 8 

System parameters for binaries with orbits 



HD 16760 HD 113449 HD 160934 



P (days) 466.5 ±0.4 216.9 ± 0.2 3764.0 ± 12.4 

To (MJD) 53336.5 ± 3 53410.5 ± 1 52389.5 ± 64 

e 0.084 ±0.003 0.300 ± 0.005 0.636 ± 0.020 

j (deg) 2.6 ± 0.5 57.5 ± 1.5 82.3 ± 0.8 

a (mas) 25.5 ± 2.8 33.7 ±0.4 152.5 ±4.7 

a (AU) 1.16 ±0.18 0.73 ±0.02 5.05 ± 0.37 

n (deg) 86.9 ±1.1 201.8 ± 1.6 266.7 ±0.6 

a; (deg) 243 ± 2 114.5 ± 0.5 216.0 ±3.1 



Mtotal (M0) 0.96 ±0.44 1.10 ±0.09 1.21 ± 0.27 



TABLE 9 
HD 160934B Astrometry 



Date 

1998 Jun 30 

2005 Apr 18 

2006 Jul 8 
2006 Sep 17 
2008 Jun 23 

2010 Apr 26 

2011 Apr 23 



MJD 

50994 
53478.9 

53924 
53995.2 
54640.3 
55312.6 
55674.6 



p (mas) 

155 ± 1 
213 ±2 
215 ±2 
218 ±2 
169.4 ±0.3 
68.8 ±0.7 
20.0 ±0.1 



S (deg) 

275.5 ±0.2 
268.5 ±0.7 
270.9 ±0.3 
271.3 ±0.7 
273.3 ±0.1 
290.0 ±0.6 
18.43 ±0.1 



Source 

Hormuth et al. (2007) 
Lafreniere et al. (2007) 
'Hormuth ct al. (2007) 
^afreniere_et_aL (2007) 

This study 

This study 

This study 



TABLE 10 
HD 160934B at fixed parallax 



7r 


Mass 


K 


Model Age 


(mas) 


(M©) 


(mag) 


(Myr) 



28.2 0.603±0.042 5.29±0.02 >100 

30.2 0.526±0.037 5.44±0.02 W^t^^ 

32.3 0.463±0.032 5.59±0.02 55±10 



TABLE 11 
Survey Detection Limits 















Detection limits for each separation annulus 
















Am (mag) 








M2 (Mj) 






Target 


Instrument 


Filter 


20-40 


40-80 


80-160 


160-240 


240-320 


20-40 


40-80 


80-160 160-240 


240-320 


AB Dor Targets 


PW And 


PHARO 


CH4s 


2.27 


4.52 


4.81 


4.74 


4.73 


381 


99 


91 


95 


96 


HIP 3589 


PHARO 


CH4s 


3.23 


5.31 


5.60 


5.54 


5.57 


443 


125 


105 


109 


107 


HIP 5191 


NIRC2 


Hcont 


4.56 


4.84 


4.79 


4.73 


4.23 


111 


103 


96 


100 


136 


HIP 6276 


NIRC2 


Kcont 


4.07 


5.33 


5.25 


5.22 


5.21 


141 


72 


75 


76 


77 


HIP 12635 


NIRC2 


Kcont 


3.85 


4.93 


4.88 


4.81 


4.83 


127 


73 


75 


78 


77 


HD 16760 


NIRC2 


Kcont 


5.16 


6.20 


6.17 


6.13 


6.13 


85 


49 


50 


50 


50 


HIP 13027 


PHARO 


Ks 


1.53 


4.38 


5.30 


5.42 


5.39 


606 


155 


96 


89 


91 


HD 19668 


PHARO 


Ks 


0.960 


4.05 


4.99 


5.09 


5.03 


609 


146 


89 


84 


87 


HIP 14807 


NIRC2 


CO 


3.98 


5.20 


5.12 


5.09 


5.08 


132 


70 


73 


74 


74 


HIP 14809 


NIRC2 


CO 


3.84 


4.95 


4.91 


4.84 


4.80 


233 


114 


118 


123 


126 


HIP 16563A 


PHARO 


Ks 


2.02 


4.71 


5.66 


5.72 


5.67 


442 


101 


58 


56 


58 


HIP 16563B 


NIRC2 


CO 


4.12 


5.26 


5.19 


5.15 


5.13 


60 


34 


36 


36 


37 


HIP 17695 


PHARO 


Ks 


0.52 


3.71 


4.74 


4.77 


4.78 


275 


57 


34 


34 


34 


HIP 18859 


PHARO 


Ks 


2.10 


4.76 


5.73 


5.80 


5.77 


632 


180 


96 


104 


94 


HIP 19183 


PHARO 


Ks 


0.32 


3.49 


4.65 


4.68 


4.69 


955 


384 


186 


183 


181 


BD+20 1790 


PHARO 


Ks 


0.89 


4.02 


5.04 


5.08 


5.09 


532 


106 


64 


63 


63 


HD 89744 


PHARO 


CH4s 


3.63 


5.69 


5.98 


5.94 


5.95 


627 


269 


223 


228 


227 


HIP 51317 


PHARO 


Ks 


2.52 


5.07 


6.01 


6.08 


6.07 


100 


28 


21 


17 


20 


HD 92945 


NIRC2 


Kcont 


4.09 


5.22 


5.15 


5.13 


5.10 


127 


71 


73 


74 


75 


GJ466 


PHARO 


Ks 


1.66 


4.44 


5.39 


5.45 


5.46 


398 


85 


51 


49 


49 


HD 113449 


PHARO 


CH4s 


2.88 


5.12 


5.52 


5.47 


5.45 


325 


90 


71 


73 


74 


EK Dra 


PHARO 


Ks 


0.80 


3.98 


4.94 


4.98 


4.98 


776 


224 


121 


117 


118 


HIP 81084 


PHARO 


CH4s 


2.23 


4.54 


4.94 


4.90 


4.90 


200 


58 


49 


49 


49 


HIP 82688 


PHARO 


CH4s 


2.71 


4.90 


5.21 


5.17 


5.16 


520 


161 


132 


136 


137 


HD 160934 


PHARO 


H 


2.11 


4.22 


4.76 


4.81 


4.64 


309 


96 


70 


68 


75 


HIP 106231 


PHARO 


CH4s 


3.18 


5.26 


5.55 


5.51 


5.46 


190 


61 


53 


54 


55 


HIP 110526 


PHARO 


CH4s 


0.67 


3.76 


3.85 


3.85 


3.72 


282 


61 


58 


58 


62 


HIP 113579 


NIRC2 


Hcont 


5.24 


5.51 


5.42 


5.36 


5.14 


98 


90 


96 


99 


104 


HIP 114066 


PHARO 


CH4s 


1.82 


4.19 


4.51 


4.50 


4.41 


271 


74 


63 


63 


66 


HIP 115162 


NIRC2 


Kcont 


4.34 


5.60 


5.51 


5.47 


5.49 


130 


67 


71 


72 


71 


HIP 118008 


NIRC2 


Kcont 


3.88 


5.11 


5.02 


5.01 


4.97 


129 


67 


71 


72 


73 


/3 Pic Targets 


























HR9 


NIRC2 


Hcont 


4.25 


4.54 


4.49 


4.41 


4.04 


134 


110 


114 


120 


154 


HIP 10680 


NIRC2 


Kcont 


3.80 


5.15 


5.05 


5.00 


4.97 


86 


30 


33 


34 


35 


HIP11437B 


NIRC2 


Kcont 


3.71 


4.91 


4.79 


4.78 


4.79 


23 


16 


16 


16 


16 


HIP 11437A 


NIRC2 


Kcont 


4.27 


5.42 


5.30 


5.30 


5.31 


27 


17 


18 


18 


18 


HIP 12545 


NIRC2 


Kcont 


3.50 


4.89 


4.84 


4.78 


4.81 


68 


21 


22 


22 


22 


51 Eri 


NIRC2 


Kcont 


5.19 


6.29 


6.21 


6.21 


6.18 


73 


24 


26 


26 


27 


HIP 25486 


PHARO 


Ks 


1.47 


4.33 


5.32 


5.37 


5.37 


489 


75 


35 


33 


33 


GJ 803 


NIRC2 


Hcont 


5.46 


5.75 


5.68 


5.64 


5.36 


16 


15 


15 


15 


16 


BD-17 6128 


PHARO 


CH4s 


3.20 


5.57 


5.57 


5.56 


5.49 


90 


20 


20 


20 


20 


HIP 112312A 


NIRC2 


Kp 


4.29 


5.31 


5.19 


5.08 


5.04 


18 


13 


13 


14 


14 


HIP 112312B 


NIRC2 


Hcont 


3.34 


3.70 


3.57 


3.53 


3.25 


19 


17 


18 


18 


20 


Her-Lyr Targets 


3 
























HD166 


PHARO 


CH4s 


2.74 


4.88 


5.13 


5.10 


5.08 


396 


126 


110 


112 


113 


HD 10008 


NIRC2 


Kcont 


3.74 


4.89 


4.80 


4.79 


4.78 


213 


112 


118 


118 


119 


HD 233153 


PHARO 


Ks 


0.54 


3.73 


4.73 


4.80 


4.80 


430 


94 


60 


58 


58 


HIP 37288 


NIRC2 


Kcont 


4.74 


5.86 


5.82 


5.78 


5.75 


66 


40 


41 


42 


43 



TABLE 11 

Survey Detection Limits 



HD 70573 


PHARO 


Ks 


0.43 


3.60 


4.58 


4.70 


4.66 


688 


227 


132 


123 


126 


HIP 53020 


NIRC2 


Kcont 


4.16 


5.34 


5.26 


5.21 


5.21 


34 


20 


20 


20 


20 


HNPeg 


PHARO 


CH4s 


2.76 


4.92 


5.17 


5.12 


5.12 


497 


174 


152 


156 


156 


Tuc-Hor Target! 


























HIP 9141 


PHARO 


Ks 


0.66 


3.86 


4.84 


4.90 


4.88 


699 


113 


60 


57 


58 


HIP 30030 


NIRC2 


Kcont 


4.64 


5.77 


5.71 


5.67 


5.63 


83 


41 


43 


44 


44 


TWA Targets 


























TWA-21 


CONICA 


L' 


1.43 


4.41 


5.38 


5.30 


5.15 


163 


19 


12 


13 


14 


TWA-6 


CONICA 


L' 


1.25 


4.34 


5.28 


5.12 


4.55 


118 


15 


8 


10 


14 


TWA-7 


CONICA 


L' 


1.64 


4.56 


5.38 


5.26 


4.02 


75 


12 


6 


7 


16 


TW Hya 


CONICA 


L' 


2.64 


5.28 


6.28 


6.20 


6.13 


79 


14 


7 


7 


7 


TWA-3 


NIRC2 


Kp 


2.43 


3.61 


3.47 


3.24 


2.96 


69 


27 


31 


40 


48 


TWA-14 


CONICA 


L' 


0.120 


3.38 


4.34 


4.23 


3.86 


316 


43 


19 


20 


26 


TWA-13B 


NIRC2 


Kp 


3.22 


4.27 


4.19 


4.03 


4.03 


73 


32 


36 


41 


40 


TWA-13A 


CONICA 


L' 


1.70 


4.57 


5.52 


5.47 


5.32 


150 


18 


12 


12 


13 


TWA-8B 


NIRC2 


Kp 


2.69 


3.89 


3.78 


3.59 


3.44 


19 


13 


14 


15 


15 


TWA-8A 


CONICA 


L' 


1.78 


4.65 


5.60 


5.50 


5.22 


81 


13 


6 


7 


9 


TWA-9 


CONICA 


L' 


1.22 


4.29 


5.32 


5.22 


4.96 


106 


15 


7 


8 


10 


TWA-23 


CONICA 


L' 


1.45 


4.44 


5.42 


5.31 


5.17 


157 


18 


11 


12 


13 


TWA-25 


CONICA 


L' 


1.56 


4.50 


5.34 


5.29 


4.79 


173 


19 


14 


14 


17 


TWA-10 


CONICA 


L' 


0.224 


3.50 


4.56 


4.43 


4.31 


190 


21 


13 


14 


15 


TWA-llB 


CONICA 


L' 


0.736 


4.05 


5.07 


5.01 


4.85 


199 


19 


12 


13 


14 


TWA-llA 


CONICA 


L' 


3.91 


6.40 


7.35 


7.23 


7.17 


160 


21 


14 


15 


15 



w 

Note. — Columns 4—8 report the 99.9% lower limits for the magnitude difference Am of any undetected companions relative to the target star obtained for each of the five separation annuli, 
which are quoted in units of milliarcseconds: namely, 20— 40 mas, 40— 80 mas, 80— 160 mas, 160-240 mas, and 240— 320 mas (see Section |6.5l >. Columns 9-13 report the corresponding upper mass P 
limits M2 for any undetected companions. The latter values were calculated using the distance and apparent magnitude of the target listed in Table[3]to first convert the magnitude differences 
in columns 4—8 into l ower limits for t he ab solute magnitude of any undetected companion. These were then converted to masses using the appropriate group age listed in Table [5] witii the 
DUSTY isochrones of lChabrier et al.l 112000') . 
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TABLE 12 
Previous Imaging 



Target Inner Angle Filter Sensitivity Minimum Mass Ref. Otlier Ref. 





(mas) 




(Amag) 


(M.i) 
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Note. — For legend to the references provided in Columns 6 and 7, see Table [3] footnotes. The quoted sensitivities are effective contrasts, 
expressed as a magnitude difference between primary and secondary, achieved by previous imaging observations in the specified filter at the quoted 
inner angle and beyond. These sensitivities are also expressed as minimum detectable companion masses over the mass range 20— 80 A/j , which 
were calculated by interpolating the DUSTY isochrones of Chabrier et al. (2000), as described in Section |6.5l However, for a subset of the targets 
we quote minimum detectable companion masses directly from the imaging paper; specifically, for the targets imaged by M05 and K03, wc quote 
the inner separation limits for detectable 10 Af j companions provided in their survey detection limit tables, and for the targets imaged by HIO, 
we quote the equivalent limits for 20 Afj companions provided in their detection limit figures. 
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TABLE 13 

fu VALUES (%) FROM MC ANALYSIS ASSUMING dM /da OC 

Full Sample Reduced Sample 

Best Upper With Her-Lyr No Her-Lyr 

djV /de QC 2e 

M power law 11 (18) 14 (24) 13 (22) 14 (22) 

g power law 10 (17) 13 (21) 12 (20) 13 (20) 

q log-normal 9 (14) 11 (18) 11 (16) 12 (17) 

M2 = 20Mj 15 (26) 21 (40) 17 (32) 19 (32) 

M2 = 40 M7 11 (19) 14 (24) 13 (22) 14 (23) 

M2 = 60Mj 9(15) 11(19) 11(17) 12(18) 

e = 

M power law 12 (17) 15 (23) 14 (21) 15 (21) 

q power law 11 (16) 14 (21) 13 (19) 14 (19) 

q log-normal 10 (13) 12 (17) 11 (16) 12 (16) 

M2 = 20Mj 16 (25) 23 (38) 18 (31) 20 (31) 

M2 = 40Mj 12 (18) 15 (24) 14 (21) 15 (22) 

M2 = 60Mj 10(14) 12(18) 11(17) 13(17) 

e = 0.9 

M power law 11 (19) 13 (24) 13 (22) 14 (23) 

q power law 10 (17) 12 (22) 12 (20) 13 (21) 

q log-normal 9 (14) 11 (19) 10 (17) 12 (18) 

M2 = 20Afj 15 (27) 20 (40) 16 (32) 18 (32) 

M2 = 40M7 11 (20) 14 (25) 13 (23) 14 (23) 

M2 = 60Mj 9(15) 11(20) 11(18) 12(19) 



Note. — Column 1 gives the values calculated using the full 67 target sample listed in Table [3] along with our "best" values for the ages and 
distances given in Tables [2] and [3] respectively. Column 2 is the same as Column 1 except that upper values for the ages and distances were used in 
the calculations, i.e. the best values plus the corresponding uncertainties listed in Tables [2] and [3] Column 3 is the same as Column 1, except that 
the 9 targets of less certain moving group membership (HD 89744, HD 92945, GJ466, EKDra, HIP 30030, TWA-21, TWA-6, TWA-14, TWA-23) 
were not included in the calculations (see Section[3j. Column 4 is the same as Column 3, except that the Her-Lyr targets are also excluded from 
the calculations. In all columns, values in parentheses are those values obtained using the aperture masking limits alone, while the values without 
parentheses are those values obtained when the detection limits from previous imaging surveys were incorporated into the calculations. 



